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ABSTRACT 


The open loop F-18 longitudinal control system is stabilized 
using H2 and H(o singular value loop shaping for a multivariable 
feedback control system. The H2 and control theories involve 
suppressing the sensitivity matrix transfer function at the lower 
frequencies for high gain performance and suppressing the 
transmissivity at higher frequencies, i.e. loop shaping. The 
singular value Bode pJot is used for MIMO systems in analogy with 
the classical Bode frequency analysis for SISO systems. There are 
two control inputs with input 1 controlling the stabilator and 
input 2 controlling the leading edge flap and trailing edge flap in 
tandem. There are two outputs: angle of attack and pitch rate. 
The Ha, design achieved a separation in that input l controlled 
angle of attack and input 2 controlled pitch rate. The first 
design is an optimum design which imposed no limitations on control 
input. A cost penalty associated with control actuator limitations 
is imposed to achieve a limited performance design. 
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I. INTRODUCTION 


Classical control analysis for single-input single-output 
(SISO) systems have the use of Bode plots, root locus technigues, 
Nyquist diagrams and simple time response analysis to judge system 
performance and stability margins. Stability margin rates the 
system's ability to withstand disturbances and/or modeling error of 
a given magnitude and still remain stable. The bandwidth in a Bode 
plot, defined as the maximum frequency at which the system response 
does not fall more than 3 db from steady state gain, is easily 
shown for SISO systems. Classical techniques usually aren't 
applicable to determining the stability margins and performance 
characteristics of MIMO systems. 

The stability margin and system performance of MIMO systems 
have been successfully evaluated using the techniques of singular 
value Bode plots of return difference matrices and loop gain 
matrices in frequency domain analysis by Doyle, Stein and Safonov 
[Refs. 1, 23 . A MIMO system is said to have good robustness if the 
system has a large stability margin, good disturbance attenuation 
and low sensitivity [Ref. 2]. 

H-infinity (Ho,) and frequency-weighted linear quadratic 
gaussian (H2 LQG) apply singular value loop shaping to MIMO 
systems. Singular value loop shaping involves shaping the system 
feedback gains over a specified frequency range in Order to meet 
system gain requirements at the lower frequencies and disturbance 








attenuation specifications at the higher frequencies. Gordon 
[Ref.3] developed numerical optimization techniques for singular 
value loop shaping which manipulates the system feedback gains as 
design parameters. 

Textbook examples of H2 and Hjo theories have been shown by 
Postrethwaite [Ref. 4] and Chiang [Ref. 5]. Chiang discusses the 
design of a hypothetical fighter design at Mach 0.9 and 25,000 feet 
altitude. Rogers and Hsu [Refs. 6 and 7] developed Ho, compensated 
designs for the X-29 for a two- and three- input longitudinal 
control system. 

The purpose of this thesis was to develop a Ho, longitudinal 
controller for the F-18 for one flight condition at Mach 0.6 at 
10,000 feet altitude. The open-loop F-18 longitudinal model was 
developed by Rojek [Ref. 8] and simplified by the author. Chapter 
II discusses the properties of MIMO feedback control systems, use 
of the return difference matrix, and the basic concepts of singular 
value loop shaping. Chapter III presents the background and 
concepts of the H2 and H,,, theories. Chapter IV describes the F-i8 
controller design beginning with the open loop F-18 uncompensated 
model, the design specifications, the design approach and the 
design results. The conclusions are given in Chapter V. The 
appendices contain the MATLAB computer code and F-18 state space 
models. 
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H. PROPERTIES OP MULTIPLE-INPUT MULTIPLE-OUTPUT 
(MIMO) FEEDBACK CONTROL SYSTEMS 

The robustness of a MIMO system includes good stability 
margin, low sensitivity to plant and controller variations, and 
good disturbance rejection to high frequency disturbance inputs. 
The above robustness properties can only be modified by altering 
the feedback paths and the associated gains. The following chapter 
will illustrate the concepts of feedback manipulation to achieve a 
good robust design upon which the H2 and H^, methods in Chapter III 
are based. Rogers [Ref. 6] and Hsu [Ref. 7] describe in detail the 
feedback properties of multivariate systems so the current chapter 
will not attempt to examine this subject in great depth. Matrices 
will be denoted by bold upper-case letters and vectors by bold 
lower-case letters in this text. 

A. MATRIX NORMS, SINGULAR VALUES AND THEIR APPLICATION TO FEEDBACK 
PROPERTIES OF MIMO SYSTEMS. 

Singular values of matrices have been found in the last decade 
to be extremely useful in extending the frequency domain Bode 
analysis of classical SISO theory to singular value Bode plots for 
MIMO systems [Refs. 1,2 and 5]. The singular values of a matrix 
A of rank r where A e c"**" are denoted by and are defined as the 
non-negative eigenvalues of A®A where ** denotes the complex 
conjugate transpose of A. The singular values are ordered such 
that a 2 >U 2 ^* • maximum singular value can be expressed 
in terms of the spectral norm: 







( 2 - 1 ) 


lA|2=max (A«A) (A) =0^ 

i 


where is the ith eigenvalue of a“a. The singular values of a 
complex nXn matrix A, a^, are the non-negative square roots of the 
eigenvalues of A**A: 

ai(A)=XY^(A»A) 

The maximum singular value 0 ^^ is given by and the minimum 
singular value is equal to a,, since the 0 's are ordered from 
aI in monotonically descending order down to 0 „. 

There are 12 useful properties of singular values listed by 
Chiang [Ref. 5) blit the three most important for the purposes of 
this paper are: 


— lAxI 

1. o (A) =roaoc xec "- - = o, 

M2 


2.fl,CA)=mih xec 


|AxI: 


M, 


^ = o, 









Property 1 is important because it establishes the greatest 
singular value of a matrix A as the maximum gain of the matrix over 
all possible directions of x. Property 2 is important because the 
least singular value of a matrix A is the minimum gain of a matrix 
A over all possible values of x. Property 3 simply states that the 
absolute value of all eigenvalues are bounded by the maximum and 
minimum singular values. 

Singular values are useful to define the maximum and minimum 
gains of the return difference matrix (to be discussed in the next 
section). For a given plant G(s) the H 2 -norm and the Hoj-norm are 
defined in terms of singular values by: 





(2-3) 


PIL A supremum a(G(j«)) 
(0 


(2-4) 


The supremum stands for least upper bound. Minimizing these norms 
form the basis of H 2 and theory. The band of maximum and 
minimum singular values plotted as a function of the frequency in 
a singular value Bode plot shovjs the degree of disturbance 
rejection, stability, and performance or system gain as reflected 
in the system bandwidth. The need to suppress high frequency plant 
disturbances suggests shaping the loop singular value so as to have 
low system gain at high frequency. The need to have a small 
sensitivity to measurement noise requires suppressing the response 
at the lower frequencies. Kwaakernak (Ref. 9) discusses for SISO 
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systems the shaping of the feedback response by reducing the scalv.r 
sensitivity S(s) at the lower frequencies and the transmissivity 
T(s) at the higher frequencies. Note that the transmissivity is 
also called the complementary sensitivity. 

B. SINGULAR VALUE LOOP SHAPING FOR MIMO ROBUSTNESS. 

Consideration of Figure 2.1 will assist in developing the 
concept of singular value loop shaping. There are two sources of 
disturbances: plant disturbance and reference noise. The 
controller is P(s) , the plant is G(s), the input is r, the control 
output from P(s) is u, and the output of the system is y. The 
output return difference matrix is I + G'-OFCs) while the input 
return difference matrix is I + F(s)G(s). The quantities F(s)G(s) 
and G(s)F(s) are the input and output loop gain matrices, 


commmd 

r — 




error 

e 


control: 


plant 

disturbance 

effects 


F(s) 


u 


I 


G(j) 




''controller” 


"plant" 


output 

_y 


Figure 2-1 MIMO Feedback Control System 
respectively. A large loop gain, L(s)=G(s)F(s) , will suppress the 
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plant disturbance but will tend to amplify measurement noise. The 
transfer matrix from d to the output y is denoted by S(s) , the 
transfer matrix from r to y is called T(s) and the transfer matrix 
from r to u is denoted by R(s). The sensitivity matrix S(s) is 
given by 

S(s) A [ I+L{s) ]-i (2-5) 

The complementary sensitivity matrix T(s) is found by 

T(s) L Lis) Il+L{s)]-^ = [I-Sis)} (2-6) 

There is no common name for R(s) which plays a part in penalizing 
the control deflections and is given as follows: 

R(s) A Fis) [J+J[/(s)]-i (2-7) 

The concept of singular value loop shaping requires high loop 
gain at low frequency to drive the sensitivity matrix S(s) to small 
values thereby suppressing plant disturbances. The sensitivity 
matrix will approach I at the higher frequencies. The 
complementary sensitivity matrix T(s) will go towards zero at 
higher frequencies. 

Since 8(s) is the closed-loop transfer matrix from the plant 
disturbance d to the plant output y, the singular values of S(s) 
determine the degree of plant disturbance attenuation of the 
system. The disturbance rejection specification is usually written 
as: 
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o (S (jw)) (jo>)) 


(2-8) 


where I is the disturbance attenuation factor. The 
disturbance attenuation factor is made a function of frequency so 
that a differeint attenuation factor can be specified for each 
frequency. 

Consideration of figure 2.2 will lead to development of the 
stability criteria as applied to MIMO systems. The effect of 


PERTURBEO PLANT 



Figure 2-2 Additive/Multiplicative Perturbations 

additive and multiplicative perturbations can be determined by the 
use of singular value plots of R{s) and T(s). The following 
discussion assumes that the system without the perturbations is 
stable. 

Jf is set to zero and taken to be the 
definition of the size of Ajj(j<i)) then the smallest Ay(s) for which 
the system is unstable is [Ref. 5] 
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(2-9) 




a(T(jo>)) 

multiplicative perturbation that the system can stand without going 
unstable. Hence the stability margin will be greater. If Aj^=0 
then the smallest A;^( jo) that makes the system unstable is 


o(A^-(ja))) = _ (2-10) 

o (R( 3 <i)) ) 

The last two equations express Robustness Theorems 1 and 2 
respectively [Ref. 5]. As a consequence of these two theorems, 
specifications can be made on the maximum allowable singular values 
of R and T matrices in terms of the weighting matrices W2(3<^^ 
W 3 (jo) as follows: 

a(R(jw)) ^ jWj'^ (j<i))| (2-11) 


a(T(j«) ) i |W3^(jo>)| (2-12) 

Figure 2-3 shows the performance boundary on a singular 
value Bode plot as determined by I Wj^(jo)l and the robustness 
boundary as determined by I W3(j o)l . The two dashed lines 
represent the minimum and maximum singular values of the loop gain 
L(s) which is the product F(s)6(s). The plot of the reciprocal of 
the maximum singular value of S(s) is seen to follow the minimum 
singular value of L(s) above the 0 db line then approaches the 0 db 
line at the higher frequencies where S(s) approaches the I matrix. 
The maximum singular value of T(s) approaches the maximum singular 
value of the loop gain L(s) below the 0 db line while approaching 
the 0 db line at the lower frequencies where T(s) approaches I. 


9 


Figure 2-3 Singular Value Specifications on 8{s) and T(s) 



III. H„ CONTROL DESIGN 
A. THEORETICAL APPROACH. 

The Hoo control problem consists of solving the small gain 
problem for the controller F(s) such that the infinity norm of the 
closed loop transfer matrix is less than or equal to 1 and is 
stable, see figure 3-1. The closed loop transfer 




yi 


Figure 3-1 The Small Gain Problem 

matrix for the plant P(s) is modified by the feedback of 
output 72 to the controller P(s) and recurning a control U 2 to the 
plant to achieve robust stability. The transfer matrix is 
defined in terms of the weighting matrices W 3 , S(s), and T(s) 
as defined by: 
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(3-1) 


T 




W^S 

W 3 T 


If one includes the weighting matrix Wj to penalize the control 
then equation (3-1) becomes: 



) 


(3-2) 


The open loop plant 6(s) must be augmented with the weighting 
matrices to form the augmented plant P(s) as shown in figure 3-2. 

The relation between the input and the output is given by 
equation (3-3). 



(3^3) 


The matrices in the brackets which contain the weighting matrices 
is the augmented plant P(s) . The state space representation of the 
augmented plant is 


/ 

A I 

Bi B 2 


'Pll I 

P12I 





_ 

CxI 

IC2! 

Dll ^12 

D 21 D22^ 


P21 ! 

i ^22, 


(3-4) 


In the H2 analysis which precedes the final design the 
matrix must be a null matrix and the Dj ^2 matrix must be full rank. 
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Augmented Plant P(s} 





Figure 3 -2 Compensated System with Augmented Plant P(s) 


The small gain problem is to find a stabilizing controller 
F(s) for the augmented plane P(s) so that the control Ujts) = 
F(s)y 2 (s) will minimize the norm of the closed loop transfer 
matrix: 

Tylul == Ph(!?)+Pi2(s) (I-I'(s)P 22 (s))'^F(s)P 2 i(s) (3-5) 

where th? closed loop transfer matrix (equation 3-1) is represented 
in terms of the partitioned matrices given in equation 3-4. 

B. DETERMINATION OF WEIGHTING CONSTRAINTS. 

The weighting constraints chosen were the same as used in the 
X-29 H„ design [Ref. 6]. The first objective is to suppress the 
sensitivity matrix singular values as much as possible for the 
.largest possible bandwidth by large loop gains. Secondly, the 
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complementary sensitivity matrix singular values must be suppressed 
by 20 db at a frequency of 100 rad/sec with a second order fall off 
of -40 db/decade for frequencies above 100 rad/sec. The resultant 
weighting constraints are: 


(Y*W,)-i(s)=Y-i* . 01=» *12X2 

(3-6) 

W2(S) = - .001*12X2 

(3-7) 

W3-^(s) = 

(3-8) 


s 


The parameter y in equation (3-6) is used in an iteration scheme 
discussed in the next section which allows us to approach the 
maximum design limits of the H2 and approaches. 




































The 0 db crossover frequency of the plot must be 

sufficiently below the 0 db crossover frequency of the W 3 plot or 
conditions (2-10) through (2-12) won't be satisfied. The V^(s) 
weighting: matrix can be seen to not have a proper state space 
representation as there are two zeros and no poles. However, 
W 3 (s)G(s) equation (3-3) for P(s) has a proper state space 
representation. The W 2 (s) weighting matrix ensures that the Djj 
submatrix has the full column rank required by theory. The 
weighting matrix was modified for the H2 design by eliminating the 
100 rad/sec corner frequency and making the denominator equal to 1 . 
The plot of the F-18 Hq, design specifications can be seen in figure 
3-3V The F-18 H2 design specifications is shown in figure 3-4. 



Frequency - Rod/Scc 

Figure 3-4 H2 Design Specifications 
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C. DISCUSSION OF H2 AND DESIGN ITERATION. 

A thorough discussion of the Hoo and H2 theory can be found in 
references 6 and 7. The H2 and Ha, techniques are usually used 
together with the H2 theory used as a first approximation with 
The Y parameter in the performance weighting matrix y^^'^Cs) is 
iteratively increased until the H2 design reaches the design 
specification limit. The final y for the H2 design is used as the 
starting point for the Hg, design whereupon the parameter y cau be 
further increased before exceeding the design constraints. This 
iterative procedure is shown in figure 3-5. 

The fact that a higher y is realizable with the Ha, design is 
an indication of the larger bandwidth, greater disturbance and 
uncertainty attenuation within the design constraints of the Hg, 
design over the H2 design. 
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Figure 3-5 H2/Ha, Design With y-Iteration 
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IV. F-18 CO:iTROL DESIGN 

A. F-^18 OPEN LOOP PITCH AXIS STATE SPACE MODEL 

The F-18 state space model is based on the model by Rojek 
[Ref. 8] for a flight condition of 0.6 Mach number at 10000 feet 
altitude. Figure 4-1 shows the F-18 with the control surfaces and 



Figure 4-1 F-18 Control Surfaces and Sign Conventions 
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the positive sign conventions for the deflections. The pitch 
control surface deflections are denoted die, dte, and dst for wing 
leading edge, wing trailing edge, and stabilator deflections, 
respectively. In the discussion of the state space representation 
to follow, the deflections will be referred to as ^st 

for the leading edge, trailing edge, p.nd stabilator, respectively. 
The multiple control inputs make the F-18 an ideal candidate for 
the robust control theory. Control input 1 controlled the 
stabilator while control input 2 controlled the leading and 
trailing flaps together. The pitch rate and angle of attack were 
chosen as the two outputs of interest so that the control system is 
a two input-two output system similar to the example of an advanced 
fighter Ha, controller treated by Chiang [Ref. 5]. 

The airframe equations of motion consists of two states: the 
downwash velocity w and the pitch rate q. Only the short period 
aircraft modes were considered neglecting the phugoid modes similar 
to Rogers [Ref. 6]. The equations of motion are linearized about 
a trim condition resulting in a set of first order differential 
equations of the general form: 

• 

X = Ax + B5 (4-1) 

Expanding the above equation in terms of the stability derivatives 
similar to that shown in McRuer [Ref. 9] gives: 
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w 

q 


^Sb -^Slf ^5tf 

1-Z 1-Z 1-Z 

M53+M^*Z53 M5if+M;,*Z5if M^tf+M^lZstf 

1-Z 1-Z 1-Z 

W 

At 10000 feet altitude and Mach 0.6 the trim angle of attack and 
corresponding pitch angle is 2.6184 degrees for level flight. Uj^ 
is the body longitudinal component of flight velocity which is 
computed knowing the true velocity of 646.42 ft/sec and the pitch 
angle. 

Figure 4-2 presents the two-input open loop actuator/ 
aircraft interface. There are two inputs U 3 ^ and U 2 with u^ being 
the input to the stabilator and U 2 the input to both the leading 
and trailing flaps. The stabilator is a fourth order actuator 
while the two flaps have second order actuators giving a total of 
eight actuator states. Scaling of the system matrix by 
transforming dbwnwash velocity w to the angle of attack by the 
relation w=V*a and transforming the units of the stabilator third 
derivative from rad/sec^ to 10^ rad/sec^ reduced the condition 
number of the system matrix from 10^ to 10^. 
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Figure 4-2 Uncompensated F-18 Open Loop Configuration 


The F-18 open loop state space model is a 10 state model 
consisting of the two airframe states and eight states for the 
three actuators. The resultant 10 state linear model of G(s)=C(si 
- A) “^B+D is presented in Appendix B. The order of the state 
variables with description and units is shown in Table 4-1 and the 
open loop poles in Table 4-2. Note there are no unstable poles in 
the F-18 open loop system matrix. The example aircraft by Chiang 
[Ref. 5] had a complex pair of unstable poles and the X-29 design 
by Rogers and Hsu [Ref. 6 and 7] had a single real unstable pole. 
The first pair of complex poles is the short period airframe poles 
with a frequency of 2.80 rad/sec. The other eight poles are the 
higher frequency actuator states. 
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4-1 The 

Ordered Uncompensated F-18 

Model States 

State 

Description 

Units 

a 

angle-of-attack 

radians 

q 

pitch rate 

rad/sec 


stabilator deflection 

rad 

is 

stabilator rate 

rad/sec 

• • 

stabilator accel. 

rad/sec^ 

• * • 

stabilator jerk 

le+04 rad/sec 


leading flap defl. 

rad 


leading flap rate 

rad/sec 

^tf 

trailing flap defl. 

rad 

^tf 

trailing flap rate 

rad/sec 


Table 4-2 

Uncompensated F-18 Open Loop Poles 

-.975 ± j 2.627 
-62.126 ± j85.022 
-^14.924 ± j33.199 
-26.902 
-82.898 

-24.850 ± j24.647 
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B. DESIGN APPROACH 

The singular values of the uncompensated F-18 open loop plant 
are plotted in figure 4-3, the upper curve is aj„ax(G(jw)) and the 



Figure 4-3 Uncompensated F-18 Open Loop Singular Value Plot 


lower curve is bandwidth of 3.7 rad/sec is narrow. 
The small loop gains for lower frequencies show 
that the F-18 open loop plant has poor disturbance rejection and is 
highly sensitive to modeling errors and system variations. From 
the H2 and Ho, control design methodologies presented, the 
sensitivity function singular values must be suppressed to the 
maximum extent possible by increasing the loop gains to as high a 















































































value for the maximum possible bandwidth without conflicting with 
the system's stability constraints. The maximum singular value 
must show an attenuation of 20 db at a frequency of 100 rad/sec 
with second order falloff of 40 db/decade. The above stability 
constraint has the purpose of attenuating the control effort at the 
higher frequencies so that the flexible structural modes ar- not 
excited. 

The weighting constraints selected for the problem are given 
in equations (3-6) through (3-8). Figures 3-3 and 3-4 showed the 
specifications for the Ho, and H2 controllers, respectively. 

The open loop plant has 10 states but the augmented plant has 
a 14th order state space representation as Wj^(s) and W 2 (s) each add 
two states to the F-18 plant 6(s) . The W 3 (s) weighting / motion, 
having no state space representation, adds no states to the 
augmented plant. The H2 and controllers will also be 14th 
order, the same as the augmented plant. 
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c. 


DESIGN RESULTS OF THE H2 AND OPTIMUM CONTROLLER 
The H2 design was undertaken first per the approach shown in 


figure 3-5 with the assumed value of 1 for y. The Matlab fl8'i2-m 
script file listed in Appendix A was used. The value of y was 
increased until the cost function IIl’yiuill 2 reached the all pass 
limit or 0 db for the H2 controller. The Hq, solution was then 
performed using the flSinf.m script file with y being increased 
until a y is reached in which any further increase will not result 
in a stabilized controller. Figures 4-4 and 4-5 show plots of the 
cost function llTyjyj|l 2 for the H2 solution with y=l and y=5.3. At 
y=l, the cT^ax singular value indicated by the solid line is about 
8 db below the all pass 0 db line. Increasing y to 5.3 pushes the 
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COST FUNCriON Tylul (Ganimii = 5.3) 



Figure 4-5 H2 Cost Function llTyiuill 2 for y=5.3 

H2 cost function l!Ty 2 yj|l 2 to the 0 db line near a frequency of 5 
rad/sec with the singular value (dotted line) pushed to within 
4 db of the 0 db line. 

In figure 4-6 the Hoo solution with a Y=i3*5 pushes the cost 
function minimum singular value to within .6 db of the 0 
db line. The significantly higher value of y possible with the Hq, 
design shows that it is clearly superior to the H2 design in 
performance. 

A comparison of the singular value plots of the sensitivity 
function S(s) and (yWj)“^(s) weighting function can be seen in 
figures 4-7 through 4-9 for the H2 design with y=l and 5.3 and the 
Hjo design with y=13.5, respectively. The singular value is the 
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Figure 4-7 Sensitivity Function S(s) for H2 Solution, y=1 
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top dotted curve in each of the figures and the singular value 

curve is represented by the dashed lower curve. In figure 4-7 the 
performance boundary levels out at -40 db below .01 rad/sec and the 
minimum singular value is 14 db below the boundary. Increasing 

Y to 5.3 in figure 4-8, the boundary dips to -54 db and the 

minimum singular value of the sensitivity function is 5 db below 
the performance boundary. The solution at A-13.5 in figure 4-9 
is now suppressed to -63 db and the sensitivity function's minimum 
singular value is now within 1 db of the performance boundary 
As Y is increased the weighting constraint (yWi)”^(s) is suppressed 
to lower magnitudes and the singular values of the sensitivity 
function 6(s) are pressed closer to the weighting function. The 
lower sensitivity curve S(s) achieved by the Hoo design shows that 
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it has superior disturbance attenuation, lower sensitivity to plant 
variations and modeling errors, and a wider control bandwidth. 

The complementary sensitivity function T(s) is plotted in 
figures 4-10 through 4-12 for the H2 case with y=1 and 5.3 and the 
Ho, design with y=13.5. The dashed curves here represent the 
singular values and the dotted curves represent the singular 

values. In figure 4-10 at y=l, the complementary sensitivity 
function, T(s), has a corner frequency of 10 rad/sec and at 100 
rad/sec the maximum singular value is at -60 db compared to the -20 
db robustness boundary. At y=5.3, the H2 solution for T{s) in 
figure 4-11 shows a corner frequency of 18 rad/sec and a gain of - 
45 db at 100 rad/sec. At y= 13*5 in figure 4-12, the solution 
yields a corner frequency of 22 db and a maximum singular value of 
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Figure 4-11 H2 Complementary Sensitivity Function T(s), y=5.3 
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Figure 4-12 Ho, Complementary Sensitivity Function, y=13.5 


-^28 db at 100 rad/sec. From these results, the higher y that the 
Hoo controller is able to achieve resulted in pushing T(s) as close 
to the robustness boundary as possible. 

The effect bn control bandwidth of increasing y can be 
quantified from figures 4-10 through 4-12. At y=l, the H2 
bandwidth determined by the frequency at the -3 db maximum singular 
value is about 10 rad/sec while the y=5.3 H2 solution achieves 
about a 16 rad/sec bandwidth. The Ho, solution at y=13.5 has a 
higher bandwidth of 20 rad/sec. Rogers [Ref. 6] found for the X-29 
that the best H2 and Hg, solutions gave 20 and 30 rad/sec, 
respectively. Both sets of results indicate that the Hq, 
compensated aircraft is a more responsive aircraft. 
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Figure 4-13 Feedback Configuration 

The controller is a 14-state controller as was expected. The 
closed loop controller as shown in figure 4-13 is a 24-state 
configuration. The output vector y consists of the output states 
a and q. The control input vector r contains the control inputs u^ 
and U 2 . Since the controller is placed in series with the F-18 
plant the commands u^ and Uj are reference commands to the outputs 
a and q. The closed loop model has 2 inputs, 2 outputs and 24 
states. 

The state space model of the 24th order closed loop Ho, 
compensated m.odel is presented in Appendix B. The poles of the Ho, 
compensated closed loop system are shown in Table 4-3. The poles 
of the open loop pla.it G(s) can be seen to also exist in the closed 
loop Ho, corpensated plant. The low condition number of the open 
loop plant G(s) means that the augmented plant P(s) and the 
controller F<s) are well-conditioned. The well-conditioned 
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Table 4-3 Ha, Compensated Closed Loop Poles 

-1.7985e+04 

-4.6415e+02 

-6.2126e+01 ± j8.5022e+01 
-1.4533e+01 ± j8.1275e+01 
-1.0183e+02 
-8.1898e+01 

-1.4924e+01 ± j3.3199e+01 
-1.9460e+01 ± j2.2795e+01 
-2.6990e+01 ± j2.3233e+01 
-2.4850e+01 ± j2.4647e+01 
-2.5633e+01 ± jl.3283e+01 
-2,6902e+01 
-3.9708e+01 

-9.7451e-01 ± j2.6269e+00 

-l.OOOOe-03 

-l.OOOOe-03 


numerical properties of these matrices required no minimum 
realization nor balancing: to be performed which is desirable as the 
meaning of the state variables becomes obscure if the matrices 
undergo balancing. The open loop 10 state matrix can be seen, in 
Appendix B, to occupy rows 15 through 24 and columns 15 through 24. 
The open loop states 1 through 10 correspond to states 15 through 
24 of the Hq, closed loop controller. 

The output return difference matrix [I+G(s)F(s)] is the 
inverse of the sensitivity matrix and as such the minimum singular 
value approximates the loop gains if the loop gains are large. The 
plot of the output return difference matrix I+G(s) of the 
uncompensated closed loop plant, figure 4-14, show that with the 
minimum singular value well below the 0 db line that the loop gain 
is low. The low loop gain of the uncompensated plant means that it 
has low disturbance rejection and high sensitivity to plant 
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Figure 4-14 Singular Value Plot l+6(s), Uncompensated F-18 



Figure 4-15 Singular Value Plot I+G(s)F(s), H^o Compensated F-18 
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variations and modeling errors as was surmised from the singular 
value plot of the open loop plant earlier in the chapter. The plot 
of the singular values of the compensated return difference matrix 
I+F(s)6(s) in figure 4-15 show that the loop gains at the lower 
frequencies are much improved over the uncompensated plant. Thus 
the Hoo compensated plant has good disturbance rejection and low 
sensitivity to plant variations and modeling errors. The steep 
second order roll-off designed into the complementary sensitivity 
function T(s) probably caused the dip of the singular values below 
0 db between 14 and 80 rad/sec indicates a lower level of 
performance near the 0 db crossover frequency. The crossover point 
of 13 rad/sec is below the 30 rad/sec Crossover of the W 3 weighting 
matrix which is one of the requirements for stability. 

The inverse-return matrix I+(G(s)F(s) )is plotted in figure 
4-16 for the Hoo design. The stability margins were determined by 
examining the universal gain and phase margin curve [Ref. 6 , pg. 68 ] 
and are the same as those guaranteed by the linear quadratic 
regulatpr problem. The minimum singular value of 0 db or 1 shown 
for the Hoo design in figure 4-16 guarantees gain margins of -6 db 
to infinity and phase margins of ±60°. 
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1. Hg, Input Return Difference Analysis 

The input additive and input multiplicative return 
difference matrices are plotted in figures 4-17 and 4-18, 
respectively. The input additive return difference matrix has poor 
disturbance attenuation with the minimum singular value (dashed 
curve) below -3 db in the range of frequency of 1 to 60 rad/sec. 
A (figure 4-18) translates to a 
gain margin of -2 to +2 db and a phase margin of 5°. In figure 4- 
18, the minimum value of the input multiplicative return difference 
matrix I+(F(s)G(s))"^ violates the W 3 boundary. The conclusion is 
that the design does not guarantee stability at the inputs to 
the plant G(s) since the design is based on the plant output. 
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Figure 4-18 F-18 SV Plot of I+(FG)^, Design for the F-18 










2. Precision Longitudinal Control Modes 

The Hoo solution in the two-input, two-output case ideally 
will result in control of output 1 by input 1 and output 2 by input 
2. This feature allows for multiple, independently controlled 
surface deflections. Safonov [Ref. 10]; listed the three precision 
longitudinal modes of control observed with the H<jo derived designs: 

1. Varying the vehicle vertical velocity by varying angle of 
attack while holding the pitch angle constant or keeping 
q equal to 0. 

2. Direct lift control by varying 6 while keeping a constant 
so that the velocity vector remains fixed along the 
aircraft stability axis Xg as Xg rotates. 

3. Pitch pointing by controlling 6 at a constant flight path 
angle so that the flight path angle or velocity vector 
remains fixed while Xg rotates (6=d). 

Closed loop Bode Plots of a and q responses to inputs 1 and 2 
are shown in figures 4-19 through 4-22. The response of a to input 
1, figure 4-19 shows a gain of 1 or 0 db for frequencies of up to 
6 rad/sec which is above the short period frequency of the F-i8. 
However the q response to input 1 in figure 4-20 never gets above 
-32 db which occurs at a frequency of 20 rad/sec, well beyond the 
short period frequency of the F-18. Figure 4-21 shows the q 
response to input 2 to have a gain of 1 or 0 db to a frequency of 
10 rad/sec while the transfer function a/U 2 in figure 4-22 is 
suppressed with a maximum of -36 db at 22 rad/seC. These Bode 
plots show the great separation with very little cross-coupling 
which allows input 1 to be used to control only a and input 2 to 
control only q. 
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Figure 4-19 Bode Plot of F-18 a/u^ T.F., Ho, Design 
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Figure 4-21 F-18 Bode Plot of q/Uj T.F. Design 
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Figure 4-22 F-18 Bode Plot of a/U 2 T-F. Design 
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The time responses of a and q to a 1° pulse of 1 second 
duration in input 1 is shown in figure 4-23 and to a l°/sec pulse 
of 1 second duration in input 2 is shown in figure 4-24. Figure 4- 
23 shows a fast rise time in the alpha response to U 2 of .2 sec to 
reach the commanded input but the q response barely makes a ripple 
along the 0 °/sec line when plotted on the same scale. Figure 4-24 
shows that the pitch rate response q slightly overshoots the 
commanded U 2 before settling out to the commanded value but here 
again the other output a is essentially zero due to the second 
input Uj. The rise time of a to u-^ is .1 sec and the rise time of 
q to U 2 is .088 sec. 
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The deflections of the stabilator, leading flap and trailing 
flap due to a 1° 1-second pulse at input 1 are plotted as a 
function of time in figures 4-25, 4-26, and 4-27, respectively. 
Since a negative deflection is required to produce a positive a 
response and the only control surface with a negative deflection 
during the l-second pulse in input 1 is the stabilator then it is 
concluded that it is the stabilator that controls a. This result 
would be expected from the open loop state space representation 
which is in the feedforward loop of the closed loop system in which 
input 1 is fed through the stabilator. 

The time response of the stabilator, leading flap, and 
trailing flap due to a 1 °/sec l-second pulse in input 2 is shown 
in figures 4-28, 4-29, and 4-30. The negative rectangular pulses 
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radians | C radians 



re 4-25 Stabilator Response to 1 /sec 1-sec Pulse in 



Figure 4-26 Leading Flap Response to 1° l-sec Pulse in Uj 







radians 



Figure 4-27 Trailing Flap Response to 1° l-sec Pulse in Uj^, H,, 
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radians 



Figure 4-30 Trailing Flap Response to l°/sec 1-sec Pulse in U 2 
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in and ^tf show that they are controlling the response to a 
positive pitch rate command U 2 and is not controlling q. 

The maximum actuator deflection limits and the actuator no 
load rate limits are given in Table 4-4; 


Table 4-4 F-T8 Actuator Deflection and Rate Limits 


Actuator Deflection Limits No Load Rate Limits 

Stabilator +10.5° 40°/sec 

-24 

Leading Flap +34 15 

- 3 

Trailing Flap +45 18 

- 8 


The maximum deflections observed occurred for input l with 2.2 rad 
for 6g, 4.4 rad for and 5 rad for The maximum deflections 
for input 2 were smaller with 1.38 rad for 5g, 2.8 rad for and 
3.1 rad for 5|.f. The actuator limits in table 4-4 were greatly 
exceeded for both actuator angular limits and rate limits. The Hjo 
solution developed here did not penalize the controls enough in the 
cost function and unrealistic actuator performance resulted. The 
Hoo limited performance design is presented next with a higher Wj to 
insure that the cost function places a greater weight on the 
controls in order to get a more realistic Hg, design. 
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D. DESIGN RESULTS OF THE LIMITED PERFORMANCE CONTROLLER 

The Ha, design was reworked to bring it within practical limits 
by increasing the W 2 e term to .018 and decreasing the corner 
frequency from 100 to 2.5 rad/sec in the weightin, function. 

The weighting function assignments are: 


(yWi(s) ) 


■ 01(1005+) 
Y(.4s+1) 


(4-3) 


W2(S) = -.018*12X2 (^"'^) 

W 3 (s) - ^ - * 12 X 2 (4''5) 

A plot of and W 3 ”^(s) weighting functions are shown in 

figure 4-31. 

The maximum y achievable with the above W2=.018I was found to 
be 1.58 which only pushed the singular value of the cost function 
to within 2 db of the all-pass 0 db line in figure 4-32. 
The optimum Ha, design with y= 1*58 in figure 4-33 shows much less 
disturbance rejection and more sensitivity to plant and modeling 
errors than the optimum Ha> design. The complementary sensitivity 
function for Y=i*58 in figure 4-34 shows a greatly reduced 
bandwidth of 2 rad/sec as against 20 rad/sec in the Ha, optimum 
design. The closed loop poles are listed in Table 4-5. 
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[ure 4-31 F-18 Limited Performance Design Specifications 
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Figure 4-32 Limited Perf. Cost Function llTy^y^lloo fot y=1.5^ 






















Figure 4-34 Limited Perf. Compl. Sensitivity Function, y=l.5Q 


49 












Table 4-5 


F-18 Limited Performance Closed Loop Poles 


-1.0017e+03 
-6.2126e+01 ± 
-1.7711e+01 ± 
-8.2898e+01 
-7.6597e+01 
-1.4924e+01 ± 
-2.4200e+01 ± 
-2.4850e+01 ± 
-2.0234e+01 ± 
-2.6902e+01 
-2.5668e+01 ± 
-9.7451e-01 ± 
-2.5001e+00 
-1.9977e+00 
-1.8000e-02 
-1.8000e-02 


j8.5022e+01 

j8.1027e+01 


j3.3199e+01 

j2.5758e+01 

j2.4647e+01 

j2.0813e+01 

j7.8886e-01 

j2.6269e+00 


Investigation of the output return difference matrices show 
how much the performance has degraded. The singular value 

figure 4-35 shows how much the performance has 
degraded. The singular value plot of [I+(G(j(i))P(jw) )in figure 
4-36 shows a steady state gain of 2 db or -7 to +3 db and a ±78° 
gain and phase margin, respectively. A typical modern fighter 
aircraft has a gain margin of -8 to +4 db and a phase margin of 
±35°. 
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Figure 4-35 Limited Perf. SV Plot of (I+GF), y=1-58 
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Figure 4-36 Limited Perf. SV plot of (l+(GF)~-^], y=1-58 
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1. Limited Performance Plant input Problem 

The singular value plot of the input additive return 
difference matrix in figure 4-37 shows that the stability of the 
limited performance F-i 8 at the plant inputs is further degraded 
over that of the optimum Ho, design. The plot of the input 
multiplicative return difference matrix I+(FG)~^ in figure 4-38 
shows that violates the W 3 boundary as did the optimum Hoo 
design. Stability robustnes is not guaranteed at the plant inputs 
as in the optimum Ho, design. 
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2. Limited Performance Precision Longitudinal Modes 

Examination of the closed loop bode plots of the four 
transfer functions a/u^/ q/u^, q/U 2 , and a/U 2 in figures 4-39 
through 4-42 shows that the cross coupling between Uj^ and q, and 
between U 2 and a is now not negligible with the limited performance 
design. The optimum had for all practical purposes no cross 
coupling. Figure 4-39 shows a steady state gain of -2 db for the 
a/input 1 transfer function which translates to an absolute gain of 
.794 instead of the 0 db gain noted for the optimum Hq, gain. The 
q/input 2 transfer function in figure 4-40 had a peak of -22 db 
instead of the -32 db of the optimum design showing an increased 
cross coupling. The results with input 2 are similar with the q/U 2 
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transfer function having a steady state gain slightly greater than 
-2 db in figure 4-41. Figure 4-42 shows the peak of the a/Uj bode 
plot to be -20 db at 2 rad/sec so not only is the peak higher than 
for the optimum Hoo case. Actually the cross coupling is a little 
stronger because the frequency at which both cross transfer 
functions peak is close to the short period frequency of 2.8 
rad/sec. The vertical scales in figures 4-40 and 4-42 are 
suppressed. 
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A 1° pulse of 1-sec duration was applied at input 1 and then 
at input 2 with the limited performance design as was done with 
the optimum design. A plot of the a and q time responses in 
figure 4-43 still show that input 1 still dominates in the control 
of angle of attack but the pitch rate response is not negligible. 
Similarly input 2 still dominates the pitch rate response but the 
angle of attack response is noticeable as seen in figure 4-44. The 
rise times are much longer as expected since the controls are being 
limited by a larger weighting value for W 2 . After 1 sec the angle 
of attack reaches .0113 radians of the .01745 radians commanded by 
input 1. The .01745 rad/sec commanded by input 2 caused the pitch 
rate to reach .0134 rad/sec after 1 sec. The rise time of the 
optimum performance design again took .2 sec in both cases. 
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re 4-43 Lim. Perf. a and q Resp. to 1 /s 1-s Pulse at u 
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The control deflections Sg, S^f, and limited 
performance F-18 are much smaller than for the optimum performance 
F-18 Hot design as seen in figures 4-45 through 4-47 for input 1 and 
figures 4-48 through 4-50 for input 2. The control deflections for 
input 1 were all within the acceptable limits for the F-18. The 
leading flap and trailing flap both exceeded the maximum negative 
deflections allowed for the F-18 of -.052 rad for the leading flap 
and -.14 rad for the trailing flap for input 2. The stabilator 
deflection was completely within its limits for either input. The 
main problem in utilization of the three control surfaces in the F- 
18 is that the leading and trailing flaps have most of their travel 
in the downward positive direction which from figure 4-1 can be 
seen to cause a negative pitching moment. The leading flap has a 
deflection range of -3 to +34 degrees and the trailing flap has a 
deflection range of -8 to +45 degrees. The X-29 had a good range 
of positive and negative travel on the canards and strakes making 
it a little more ideal for the Ho, approach [Ref. 6]. NASA Dryden 
has been testing a more maneuverable F-18 modified to use canards 
which would make it a more ideal candidate for the Hq, approach. 
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Figure 4-45 F-18 Lim. Perf. Stabilator Sg to 1° 1-s Pulse at Uj 
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Figure 4-47 Lim. Perf. Trailing Flap to 1° 1-s Pulse at u^ 
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Figure 4-50 Liir.. Perf. trailing Flap 5^.^ to l°/s Is Pulse at U 2 
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The actuator no load rate limits shown in table 4-4 were 
exceeded in both input 1 and input 2. The control rates of the 
actuators are shown in figures 4-51 through 4-53 for the pulse from 
input 1 and figures 4-54 through- 4-56 for the pulse from input 2. 
The peak rates range from .8 to 6 rad/sec which exceed the 
capabilities of the actuators. 
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'igure 4-52 Lim. Perf. d5j^f/dt due to 1° Is Pulse in u^ 



Figure 4-53 Lim. Perf. d5j.£/dt for 1° Is Pulse in Uj^ 
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radians/second 1 radians/second 



Figure 4-54 Lim. Perf. d5g/dt for l°/s Is Pulse at U 2 
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F-18 Uni. D'n-F:>OT FOR O.On-IS rad / 1 .-icc STEI’ (INPUT 2) \V2=--..()I8 



Figure 4-56 Lim. Perf. d5^.f/dt due to l°/s Is Pulse at U 2 

The approach to a limited control design did show a reduction 
in control deflections and rates by increasing the control weight 
in the cost function but at the cost of system performance and 
robustness. In addition to the higher control weighting the F-18 
weighting performance specification had to be relaxed. The 
small control weightings in the optimum design required 

unacceptably high energy influx into the system beyond the 
capability of the actuators to provide. The limited performance 
design still demonstrated predominantly separate control of a by 
input 1 and q by input 2 but not near the separation demonstrated 
by the optimum design. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

Hoo control theory implemented with the Matlab Robust Controls 
Toolbox was shown to be a systematic and straightforward method to 
shape the frequency response of a multi-output, multi-output system 
of the F-18 to achieve both robustness and system performance. The 
decoupling of the pitch rate and angle of attack states 
demonstrates the potential for precision flying modes in which 
different independent control surfaces control different outputs. 
The stabilator was controlled by input 1 and was able to control 
angle of attack separately from input 2 which controlled the pitch 
rate by commanding the leading and trailing flaps which were ganged 
together. It was shown that the performance level of the optimum 
performance Hq, design placed too great of an energy demand on the 
actuators. A limited performance Hoo design was demonstrated which 
placed a greater weight on the control energy in the Hq, cost 
function. The sensitivity weighting function specification had to 
be relaxed to arrive at the final design. The control deflection 
rates are still too high and the design will have to be limited 
still further. 

The Hoo theory was shown to be superior to the H2 frequency 
weighted linear quadratic regulator theory by achieving higher 
levels of plant disturbance attenuation, better suppression of 
plant variations and modeling errors, and wider system control 
bandwidth. Stability of the plant outputs was demonstrated 
although it was found that the plant inputs did not meet the design 
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specifications. The inputs failed the design specifications 
because the specifications are in tenas of the sensitivity and 
complementary sensitivity functions which are formulated in terms 
of the plant outputs. 

Actuator performance has been found to be the most severely 
limiting factor in the Ha, design and has to be more directly 
accounted for. Including the control deflections and/or rates in 
the output vector as Rogers [Ref. 6] suggested would be a way to 
make the design process more practical. 

The promise of Ha, design theory may make it more practical to 
include higher performance actuators in future aircraft design. 
Faster actuators mean cost penalties and more weight. Improving 
the speed of the actuators must be shown to improve the 
capabilities and performance of the aircraft sufficiently to pay 
for the cost and weight penalties. A concern of using faster 
actuators in the past has been the excitation of flexible 
structural modes. An Ho, design which can suppress the flexible 
structural modes, provide independently controlled outputs, 
suppress disturbances and plant variations and achieve respectable 
control bandwidths might provide the incentive to use improved 
actuators. 

Applying Ho, design to the highly maneuverable modified F-18 
being tested at NASA Dryden with canards and flow deflectors would 
seem to be a worthwhile project for future investigation. The Ho, 
design analysis should be applied to other F-18 flight conditions. 
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APPENDIX A 


MATLAB SCRIPT FILES 


68 



fl8h2.ro 

diary fl8h2.dat 
forroat short e 
disp(' ') 

disp(' ') 

disp(' This script file is designed to solve the H2 optiroal 

control ') 

dispC problem for the F-18. The 10th order FDLTI roodel, in 

state ') 

dispC space form, is that of the F-18 aircraft and actuator 

dynamics.•) 

dispC Four states are those of the aircraft dynamics, i.e., 

alpha & q') 

dispC The remaining 8 states are the dynamics of a fourth ') 

dispC order and two 2nd order actuators, i.e., the stabilator, 

leading ') 

dispC flap and trailing flap actuators. The order of the 

unbalanced ') 

dispC states is as follows;') 

dispC ') 

dispC alpha, q, ds, dsdot, dsdbldot, dstrpldot, dlf, 

dlfdot,•) 

dispC dtf, dtfdot •) 

disp C •) 

disp(' Given the open loop transfer function G(s)=Cinv(Ts-A)+D, 
a •) 

dispC stabilizing controller F(s) will be found such that the 
H2 norm') 

dispC of Tylul is minimized. ') 

disp(' ') 

disp(' H2 optimal control synthesis is performed to determine 

attainable') 

dispC performance levels. Once completed, an Hinf optimal 

control ') 

dispC synthesis is performed.') 

disp(' ') 

disp(' ') 

%pause 

clc 

disp(' The scaled F-18 aircraft and actuator state space 

representation') 
disp(' ') 

disp(' ') 

ag=[-.114d+01 9.7938d-01 -9.9461d-02 0.0 0.0 0.0 ... 


.6206d- 

■04 0. 

0 -3. 

9669d- 

-03 0 

.0;-7.0738d-H00 -.80902d-»-00 ... 

-7.9063d-00 

0.0 

0.0 

0.0 

-3.008d-02 0.0 -1.8786d-00 0.0; 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 0.0 0.0 0.0 0.0; 

0.0 

0.0 

0.0 

0.0 

1.0 

0.0 0.0 0.0 0.0 0.0; 

0.0 

0.0 

0.0 

0.0 

0.0 

ld-f04 0.0 0.0 0.0 0.0; 

0.0 

0.0 - 

•.14691d-f04 

-.49559d-^02 -.16122d-f01 -.1541d-»-03 0.0 0.0 

O 

• 

o 

• 

• 

o 

• 

o 

O 

• 

o 

o 

• 

o 

o 

• 

o 

0.0 0.0 0.0 0.0 1.0 0.0 0.0; 
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0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 -.1225d+04 -.497d+02] 


bg=[0.0 0.0; 

0.0 0 . 0 ; 

0.0 0 . 0 ; 

.21377d+04 0.0; 

-.30532d+06 0.0; 

.27277d+04 0.0; 

0.0 0 . 0 ; 

0.0 .22301d+04; 

0.0 0 . 0 ; 

0.0 .1225d+04] 

cg=[1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 

0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0] 

dg=zeros(2) 

%pause 
disp(' *) 


dispC •) 

disp(' Calculate the poles and transmission zeros of the open 
loop plant') 
disp(' ') 

poleg=eig(ag ), tzerog=tzero(ag,bg,cg,dg) 
dispC ') 

disp(' ') 

disp(' Determine the condition number of ag') 

disp(' ') 

disp(‘ •) 

condag=cond(ag), rcondag=rcond(ag) 
disp(' ') 

dispC ') 

%pause 
dispC ') 

dispC « Design Specifications » ') 

dispC ') 

dispC 1) • Robustness Spec. : -40 dB roll-off, -20 db @ 100 



-1 


-1 (100s + 1) 


disp(* Wl(s) = Gam *-* i') 

disp(' 100 2x2') 

disp(' ') 


dispC where "Gam" in this design is iteratively updated 

from i') 

w=logspace(-3,3,100); 

)c=1000; mn=[2 2]; tau=0.0; 

nuw3i = [0.0 k]; dnw3i = [1.0 0 0]; 

svw3i = bode(nuw3i,dnw3i,w); svw3i = 20*logl0(svw3i); 

nuwli = [100.0 1.0]; dnwli =[0 100.0]; 

svwli = bode(nuwli,dnwli,w); svwli = 20*logl0(svwli); 

aw2=-.001*eye(2); bw2=zeros(2); cw2=zeros(2); dw2=-.001*eye(2); 

dispC •) 

dispC •) 

disp(' (strike a key to see the plot of the 

weightings ...)•) 

%pause 

semilogx(w,svwli,w,svw3i) 
grid 

title('F-IS Design Specifications') 
xlabel('Frequency - Rad/Sec') 
ylabel('l/Wl & 1/W3 - dB') 
text(.01,0,'1/Wl(s)•) 
text(.5,100,'1/W3(s)') 
meta flSspecs 
%pause 
clc 

disp(' « Problem Formulation »') 

dispC ') 

disp(' Form an augmented plant P(s) with these two weighting 
functions:') 
dispC •) 

dispC 1) • Gam*Wl penalizing error signal "e"') 

dispC •) 

dispC 2). W3 penalizing plant output "y"') 

disp C ') 

dispC and find a stabilizing controller F(s) such that the H2 
norm') 

dispC of TF Tylul is minimized, i.e.') 
dispC ') 

dispC roin |Tylul| < i,*) 

dispC F(s) 2 ') 

dispC •) 

disp(' where ') 

dispC ! -11') 

dispC Tylul = j Gam*wi*(l + GF) | = J Gam * W1 * S 

I') 

dispC 1 -1| 1 W3 * (I - S) 

1 ') 

dispC 1 W3*GF*(I + GF) j') 

dispC •) 
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disp(' 

•) 


dispC 
disp(' 
...)•) 
%pause 
clc 

') 


dispC 
disp(' 
dispC 
disp(' 
dispC 

•) 


•) 


•) 



k 

* 

disp(' 

•k * -k * 


AUGMENT.M 

or 

disp(' 

*’) 

k 

disp(' 

*') 

k 

disp(' 


k 


*') 


disp(* 

*') 

k 

dispC 

*•) 

k 

disp(' 

*') 

k 

dispC 

*') 

k 

disp(' 

*') 

k 

disp(' 

*') 

k 

disp(» 


k 

* * * ■ 

* * * *) 



(strike a key to continue 


« DESIGN PROCEDURE »•) 
*********************** 

* [Step 1]. Do plant augmentation (run 

*.) 

AUGF18.M) 

[Step 2]. Do H2 synthesis (run H2LQG.M) 

[Step 3]. Redo the plant augmentation for a 
new "Gam" and rerun H2LQG.M 

[Step 4]. Redo the plant augmentation for a 
higher "Gam” then run HINF.M 
*********************** 


') 

•) 


disp(' 
dispC 
disp(' 

...)•) 

%pause 
clc 
disp(' 
dispC 
dispC 
dispC 

disp(' ') 

disp(' 
pack 

gama=[6.7]; 
ngaraa=length(gama); 
for i=l:ngama 
Gam=gama(l,i) 


’) 

•) 

') 


(strike a key to continue 


Assign the cost coefficients "Gam" with Gam=l ') 
serving as the baseline design 


k=1000 
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disp(' ') 

d i s p ( ' 


disp(' 

dynamics') 
d i s p ( ' 


- .) 

augment % Plant augmentation of the F-18 


- ,) 

disp(' •) 

disp(' •) 
augflS 
disp(' ') 

disp(* •) 

disp(' (strike a key to continue 

...)•) 

[ac,bc,cc,t,kl]=ctrbf(A,[B1 B2] , [C1;C2]); 

sum(kl) 

pause 

clc 

%disp(' ') 

%disp(' •) 

%disp(' Do state space balancing on the augmented plant if 

needed') 

%disp(• •) 

%disp(• •) 

%[augA,augB,augC,augD]=roinreal(A, [B1B2], [C1;C2], [DllD12;D21 D22]) 
%[abal , bbal,cbal,g,t]=obalreal(augA,augB,augC) 
%[abal,bbal,cbal,g,t]=obalreal(A,[B1 B2],[Cl;C2]) 

%A=abal, Bl=bbal(:,1:2), B2=bbal(:,3:4), Cl=cbal(l:6,:), 
C2=cbal(7:8,:) 
disp(' ') 

disp(' ') 

disp(' The transmission zeros, poles and condition number of the 
augmented') 

dispC plant follow. In addition, determine if (A,B1) & {A,B2) 
are ') 

dispC stabilizable and if (Cl,A) & (C2,A) are detectable.') 
disp(' ') 

disp(* ') 

%tzeroaug=tzero(abal,bbal,cbal,dbal) 
tzeroaug=tzero(A,[B1 B2],[C1;C2],[Dll D12;D21 D22]) 
poleaugA=eig(A) 

%poleaugA=eig(abal) 
condaugA=cond(A) 
rcohdaugA=rcohd(A) 

%condaugA=cond(abal) 

%rcondaugA=rcond(aba1) 
eps=eps 

toldef=ld*max(size(A))*norm(A,1)*eps 
%toldef=10*max(size(abal))*norm(abal,1)*eps 
tol=l00*eps*norm([A Bl]) 

[Ale,Blc,Clc,t,k]=ctrbf(A,Bl,Cl,tol) 
tol=i00*eps*norm{[A B2]) 
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[A2C,B2c,C2c,t,k]=ctrbf(A,B2 , C2,tol) 
tol=100*eps*norm([A;C1]) 

[Alo,Bio,CIO,t,k]=obsvf(A,Bl,Cl,tol) 
tol=100*eps*norm([A;C2]) 

[A2o,B2o,C2o,t,k]=obsvf(A,B2,C2,tol) 

clear condag rcondag poleg tzerog svw3i svwli 

clear condaugA rcondaugA 

clear functions 

%pause 

disp(' ') 

disp(* *) 

d i s p ( • 

----- .) 

dispC h2lgg % Running script file H2IiQG.M for H2 

optimization*) 

d i s p ( * 

------- .) 

aretype='Schur' 
h21gg 
disp(' ') 

disp(' ') 

disp(' (strike a key to continue 

...)•) 

clear functions 
%pause 

%pltopmod % Preparing singular values for plotting 

%end 

disp(' ') 

disp(» •) 

disp(' (strike a key to continue 

%pause 
disp(' ') 

dispC •) 

disp(' State space representation of controller (acp, bcp, ccp, 
dcp) •) 

dispC CLTF Tylul (acl, bcl, ccl, del) and poles, 

controllability, ') 

dispC observability, and condition number of the controller. ') 
dispC •) 

dispC ') 

pltopmod %preparing singular values for plotting 

end 

acp,bep,cep,dcp 
polec=eig(acp) 
tol=100*eps*norm((acp bep]) 

(aepe,bepe,cepe,t,k]=ctrbf(acp,bep,cep,tol) 

lc=sum(k) 

pause 

tol=100*eps*norm([acp;ccp]) 

[aepo,bepo,cepo,t,k]=obsvf(acp,bep,cep,tol) 
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lo=suni(k) 

pause 

condacp=cond(acp), rcondacp=rcond(acp) 

acl,bcl,ccl,del 

polet=eig(acl) 

disp(' ') 

disp(' ') 

disp(' (strike a key to continue 

...)') 

%pause 
disp(• ') 

disp(' ') 

disp(' Open loop state space representation of controller/plant 
series') 
disp(' ') 

[algf,blgf , clgf , uigf]=series(acp,bep , cep , dep,ag,bg,eg,dg) 
polol=eig(algf) 
disp(' ') 

disp(' ') 

disp(' (strike a key to continue 

...)') 

%pause 
disp(' •) 

disp(' ') 

disp(' Closed loop state space representation of 

controller/plant series,') 

dispC controllability, observability, and condition number of 
the closed') 

disp(' loop aegf matrix. ') 
disp(' ') 

disp(' ') 

[aegf,begf,ccgf,degf]=feedbk(algf,blgf,clgf,dlgf,2) 
tol=100*eps*norm([aegf begf]) 

[ n.cgf c, begf c, ccgf c, t, k] =ctrbf (aegf, begf, ccgf, tol) 
tol=100*eps*norm([aegf;ccgf]) 

[aegfo,begfo,ccgfo,t,k]=obsvf(aegf,begf,ccgf,to1) 
condaegf=cond(aegf) 
disp(' ') 

disp(' ') 

disp(' (strike a key to continue 

...)•) 

%pause 
disp(' ') 

disp(' Poles of the closed loop system') 
disp(• ') 

polcl=eig(acgf) 
end 
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%fl8hinf.m 
diary fl8hinf.dat 
format short e 
dispC fl8hinf,m') 
dlsp(' ’) 

disp(' This script file is designed to solve the Hinf optimal 

control *) 

disp(* problem for the f-18. The 10th order FDLTI model, in 

state ') 

dispC space form, is that of the f-18 aircraft and actuator 

dynamics.') 

dispC Two states are those of the aircraft dynamics, i.e., 

alpha & q ') 

disp(' The remaining 8 states are the dynamics of the fourth 

order•) 

dispC stabilator, the 2nd order leading flap, and the 2nd 

order trailing') 

dispC flap actuators. The order of the unbalanced states is 

as ') 

disp(' follows: ') 

disp C ') 

dispC alpha, q, ds, dsdot, dsdbldot, dstrpldot,dlf,dlfdot,' ) 

dispC dtf ,dtfdot') 

dispC ') 

dispC Given the open loop transfer function G(s)=Cinv(Is-A)+D, 
a ') 

dispC stabilizing controller F(s) will be found such that the 
Hinf norm') 

dispC of Tylul is less than or equal to one. ') 

disp(' ') 

dispC ') 

pause 

clc 

dispC F-18 aircraft and actuator state space representation') 
disp(' ') 

ag=[-.114d+01 9.7938d-01 -9.9461d-02 0.0 0.0 0.0 ... 

5.6206d-04 0.0 -3.9669d-03 0.0;-7.0738d+00 -.80902d+00 ... 


-7.9063d-00 

0.0 

0.0 

0.0 

-3.008d-02 0.0 

1 -1.8786d-00 0.0; 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 

0.0 0.0 

0.0 0.0; 

0.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 0.0 

0.0 0.0; 

0.0 

0.0 

0.0 

0.0 

0.0 

ld+04 

0.0 0.0 

0.0 0.0; 

0.0 

0.0 - 

■.14691d+04 

-.49559d+02 -.16122d+01 -.1541d+03 0.0 0.0 

• • • 

0.0 

0.0; 

0.0 

0.0 

0.0 

0.0 

0.0 0.0 

0.0 1.0 0.0 0.0; 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 - 

-.22301d+04 

-.1098d+03 0.0 0.0; 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 0.0 

0.0 1.0; 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 0.0 - 

•.1225d+04 -.497d+02 3 ; 

ag(2,7)= 

=ag(2. 

7)*57 

.296 






bg=[0.0 0.0; 

0.0 0 . 0 ; 
0.0 0 . 0 ; 
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.21377d+04 0.0; 

-.30532d+06 0.0; 

.27277d+04 0.0; 

0.0 0 . 0 ; 

0.0 .2236ld+04; 

0.0 0 . 0 ; 

0.0 .1225d+04] 

[ 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 . 0 ; 

0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 . 0 ] 

dg= 2 eros(2) 

%pause 
disp(' •) 

poleg=eig(ag), tzerog=tzero(ag,bg,cg,dg) 
dispC ') 
disp(' ') 

disp(' Determine determine the condition number of ag *' 
disp(' ') 
disp(* ') 

condag=cond(ag), rcondag=rcond(ag), rank(ag), 
pause 

disp(' ') 

disp(' ') 


%pause 
disp(' ') 
disp(• 
disp (' '). 
disp(» 1) 

Rad/Sec.') 
disp(' 

disp(’ ’) 

disp(' 
disp(' 
disp(' 
disp(' 
disp(' ') 
dispC •) 
disp(' 2) 

function') 
disp(' 


« Design Specifications » ') 

Robustness Spec. : -40 dB roll-off, -20 db @ 100 
Associated Weighting:') 

-1 1000 ') 

W3 (s) = -- * I (fixd) ') 

2 2x2') 

s ') 


Performance Spec.: minimizing the sensitivity 

as much as possible.') 
Associated Weighting:') 

') 

-1 -1 .01(100s + 1) ') 

Wl(s) = Gam *-* I') 

2x2') 

(.Ols + 1) ') 


disp(' 
disp(' ') 
disp(' 
disp(' 
disp(• 
disp(' 
disp(' 
disp(' ') 

disp(' note gamma 13.5 .Ols + 1 and limited gamma 1.58 .4s + 

1 ') 

disp(' where "Gam" in this design is iteratively updated 

from 1') 
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w=logspace(-3,3,100); 

k=1000.; inn=[2 2]; tau=0.0; 

nuw3i = [0,0 k]; dnw3i = [1.0 6.0 0.0]; 

svw3i = bode(nuw3i,dnw3i,w); svw3i = 20*logl0(svw3i); 

nuwli = [1 .01]; dnwli =[.01 1.]; 

svwli = bode(nuwli,dnwli,w); svwli = 20*logl0(svwli); 

dispC w2 weighing matrix 13.5 =-.001 1.55 = -.018') 

aw2=-0.001*eye(2); bw 2 = 2 eros( 2 ); cw2=zeros(2); dw2=-0.001*eye(2); 

dispC •) 
dispC ') 

disp(' (strike a key to see the plot of the 

weightings ...)') 

%pause 

semilogx(w,svwli,w,svw3i) 
grid 

title('F-18 Design Specifications-') 
xlabel('Frequency - rad/sec') 
ylabel('l/Wl & 1/W3 - db') 
text(.01,0,'l/Wl(s)') 
text(.5,100,'1/W3(s)') 

%meta flSspecs 
%pause 

clc 

disp(' « Problem Formulation »') 

disp(' ') 

disp(' Form an augmented plant P(s) with these two weighting 
functions:') 
disp(' ') 

disp(' 1). Gam*Wl penalizing error signal "e"') 

disp(' ') 

disp(' 2). W3 penalizing plant output "y"') 

disp(' ') 

disp(' and find a stabilizing controller F(s) such that the Hinf 
norm') 

dispC of TF Tylul is less than or equal to one, i.e.') 
dispC ') 

dispC jTylul{ < or = 1,') 

disp(' F(s) inf) 

disp(' ') 

disp(' where ') 

dispC I -IC) 

dispC Tylul = { Gam*Wl*(I + GF) | = | Gam * W1 * S 

! ') 

dispC j -l! j W3 * (I - S) 

1 ') 

dispC I W3*GF*(I + GF) j') 

disp(' ') 
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(strike a key to continue 


dispC V) 
disp(' ') 
disp(' 

%pause 
clc 

dispC ') 
disp(' ') 

dispC « DESIGN PROCEDURE»') 

dispC ') 

disp C ************************ 

* * * * * *•) 

dispC * [Step 1]. Do plant augmentation (run 

AUGMENT.M or *') 

disp(* * AUGX29.M) 

*') 

disp(• * 

*.) 

dispC * [Step 2]. Balance the augmented plant for 

better *') 

dispC * numerical condition if necessary 

*.) 

disp(' * 

*') 

disp(‘ * [Step 3]. Do Hinf synthesis with "Gam" = 1 

*•) 

disp(' * 

*.) 

dispC * [Step 4]. Redo the plant augmentation for a 

*.) 

dispC * higher "Gam" and rerun HINF.M 

* •) 

dispC ************************ 

******•) 
dispC ') 
dispC •) 

dispC (strike a key to continue 

...)•) 

%pause 

dispC •) 
dispC •) 
disp C ') 

dispC Assign the cost coefficients "Gam" with Gam=l ') 

disp C ') 

disp(• serving as the baseline design 

gama=[10.5]; 

hgaroa=length(gama); 

for i=l:ngama 

Gam=gama(l,i) 

k=1000 

dispC •) 
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d i s p ( ' 


--- •) 

dispC augment % Plant augmentation of the F-18 

dynamics•) 
d i s p ( ' 

----- .) 

dispC •) 
disp(* ') 

augflS; 
disp(' ') 

disp(' ') 

disp(' (strike a key to continue 

...)') 

%pause 

dispC ') 
disp(' *) 

%dispC Do state space balancing on the augmented plant if 

needed•) 

%disp(' •) 

%disp(' ') 

%[abal,bbal,cbal,g,t]=obalreal(A,[B1 B2],[C1;C2]) 

%A=abal, Bl=bbal(:,1:2), B2=bbal(:,3:4), Cl=cbal(l:4,:), 
C2=cbal(5:6,:) 
disp(' ') 
dispC ') 

dispC The transmission zeros, poles and condition number of the 
augmented') 

disp(* plant follow. In addition, determine if (A,B1) & (A,B2) 
are ') 

dispC stabilizable and if (Cl,A) & (C2,A) are detectable.') 
dispC ') 
disp C ') 

tzeroaug=tzero(A,[B1 B2],[C1;C2],[Dll D12;D21 D22]), 

poleaugA=eig(A) 

condaugA=cond(A), rcondaugA=rcond(A) 


eps=eps 

toldef=10*max(size(A))*norra(A,1)*eps 
tol=100*eps*norm([A Bl]) 

[Alc,Blc,Clc,t,k]=ctrbf(A,Bl,Cl,tol) 
tol=100*eps*norm([A B2]) 

[A2C,B2c,C2c,t,k]=ctrbf(A,B2,C2,tol) 
tol=100*eps*norm([A;C13) 

[Alo,Bio,CIO,t,k]=obsvf(A,Bl,Cl,tol) 

tol=100*eps*norm([A;C2]) 

(A2o,B2o,C2o,t,k]=obsvf(A,B2,C2,tol) 

clear condag rcondag poleg tzerog svw3i svwli 

clear condaugA rcondaugA 

clear functions 

%pause 
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disp(* ') 

disp(' ') 

d i s p ( ' 

----- .) 

dispC hinf % Running script file HINF.M for Hinf 

optimization') 

d i s p ( • 



-- .) 

aretype='Schur' 
hinf 

disp(' ') 

disp(• ') 

disp(' (strike a key to continue 

...)') 

clear functions 
%pause 
disp(• •) 

disp(• •) 

disp(' State space representation of the full order controller*) 
disp(' (acp, bcp, cep, dep) with its poles and condition 

number') 
disp(• •) 

acp,bep,cep,dep 
polec=eig(acp) 
condacp=cond(acp) 
rank(acp) 
pause 

toldef=10*max(size(acp))*norm(acp,1)*eps 
tol=100*eps*norm([acp bcp;ccp dep]). 

disp(' Poles, controllability, observability, and condition of 

the ') 

disp(' controller') 

disp(' ') 


polet=eig(acl) 

%pause 
disp(' 
disp(* 
pltopmod 
end 
dispC 
disp(• 
disp(* 

...)') 

%pause 
disp(* 
disp(' 

disp(' Open loop state space representation of controller/plant 
series') 
disp(' ') 

[algf,blgf,clgf,dlgf]=series(acp,bep,cep,dep,ag,bg,eg,dg) 


') 

') 

•) 

•) 


') 

’) 


% Preparing singular values for plotting 


(strike a key to continue 
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polol=eig(algf) 
dispC ') 
disp(' ') 

disp(' (strike a key to continue 

...)•) 

%pause 
dispC ') 
disp(' ') 

disp(• Closed loop state space representation of 

controller/plant series,*) 

disp(' controllability, observability, and condition number of 
the closed') 

dispC loop acgf matrix. *) 
disp(' ') 

disp C ') 

[acgf,bcgf,ccgf,dcgf]=feedbk(algf,blgf,clgf,dlgf,2) 
tol=lOO*eps*norm([acgf bcgf]) 

[acgfc,bcgfc,ccgfc,t,k]=ctrbf(acgf,bcgf,ccgf,tol) 
tol=100*eps*norm([acgf;ccgf j) 

[acgfo,bcgfo,ccgfo,t,k]=obsvf(acgf,bcgf,ccgf,tol) 
condacgf=cond(acgf) 
disp(* ') 

disp(' •) 

disp(' (strike a key to continue 

...)•) 

%pause 
dispC ') 

disp(' Poles of the closed loop system *) 
dispC •) 
polcl=eig(acgf) 
end 
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%augf18.m 

% Plant Augmentation for the F-18 H2 and Hinf problem as W3 is not 
a 

% proper transfer function. Includes contingency for adding W2 to 
% ensure dl2 is full column rank. This script file is designed for 
% the F-18 system without theta as a state, ie, 2 inputs & 2 
outputs. 
disp(' •) 

disp(' ') 

disp(' « Plant Augmentation »•) 

%Gam=gama(1,i) 

Gam = input(' Input the cost coefficient "Gam'' = '); 

cgb = l/k*[cg(l,:)*ag*ag;cg(2,:)*ag*ag*ag*tau+cg(2,:)*ag*ag] 
dgb = l/k*[cg(l,:)*ag*bg;cg(2,:)*ag*ag*bg*tau] 
nwl = Gam*[dnwli;0 0;0 0;dnwli] 
dwl = nuwli 

sysw2=[aw2 bw2;cw2 dw2]; xw2=2; 

%sysw2=[]; xw2=0; 

%[A,Bl,B2,Cl,C2,Dll,D12,D21,D22]=augmod(ag,bg,cg,cgb,dg,dgb,nwl,d 

wl,mn) 

[A,B1,B2,C1,C2,D11,D12,D21,D22]=augf18pl(ag,bg,eg,cgb , dg , dgb,nwl , 
dwl,sysw2 , xw2,mn) 

disp(' ') 

disp(' _ « _ State-Space (A,Bl,B2,Cl,C2,Dll,D12,D21,D22) is 

ready for') 

dispC the Small-Gain problem - - -•) 

% 

%-End of AUGF18.M-% 


83 













function 
[a,bl,b2,cl,c2/dll,dl2,d21,d22]=augf18pl(ag,bg,cga/cgb , dga,dgb,nu 
in, den / sy s, X, inn) 

% 

% 

[A,B1,B2,C1,C2,D11,D12,D21,D22]=AUGF18PL(AG,BG,CG,CGB,DG,DGB,NUM, 

DEN,SYS,X,MN) 

% augf 18pl produces the augmented plant for the F-18 h2 and hinf 

problem for 

% an improper W3 weighting and W2 added to make dl2 full 

column rank. 

% 

% Weighted Plant : 

% 

% Xg = Ag Xg + Bg Ug 

% 

% jYgaj = ICuc-j Xg + !Dga| Ug 

% lYgbl [Cgt: iDgbl 

% 


% 

Weighting Wl : 



Weighting W2 



% 

• 




• 




% 

% 

% 

% 

% 

Xwl = 

Awl Xwl + 

Bwl Uwl 

Xw2 

= Aw2 XW2 + Bw2 

Uw2 

YWI = 

Cwl Xwl + 

Dwl 

Uwl 

yw2 

= Cw2 

Xw2 + Dw2 

Uw2 

Over all 

augmented plant ; 





% 

• 








% 

|xg 1 

1 Ag 

0 

0 1 

ixg i 

!o 

Bg j 

lui! 

% 

1 1 _ 

1 ♦ , — 

1 

1 


1 

ixwij 

+ 1 

1 

1 

iU2i 

% 

% 

% 

% 

% 

ixwij 

1 1 

1-BwlCga 

1 

Awl 

0 ! 

1 

|Xw2j 

jBwl 

1 

-BwlDga1 

1 


1 * 1 
|Xw2i 

1 

! 0 

0 

AW2 1 


1 

!o 

1 

Bw2 i 


lYlaj 

1-DwlCga 

Cwl 

0 ! 

!xg 1 

jowl 

-DwlDga1 

lull 

% 

! Ylbj 

! 0 

0 

CW2 j 

|xwij 

jo 

Dw2 1 

1U21 

% 

lYlcj = 

1 Cgb 

0 

0 I 

|Xw2j 

+ |o 

Dgb j 


% 

!Y2 1 

1 -Cga 

0 

0 ! 


11 

-Dga 1 


% state 

Space of Weighting 

Wl & 

W2: 






% 

[awl^bwl,cwl,dwl] = tfm2ss(num,den,mn(l),mn(2)) 

[aw2,bw2,cw2,dw2] = sys2ss(sys,x) 

% 

% State Space of Augmented Plant : 

% 

%-A matrix : 

% 

[rag^cag] = size(ag); 

[rawl,cawl] = size(awl); 

[raw2,caw2] = size(aw2); 

a = [ag zeros(rag,cawl) zeros(rag,caw2);-bwl*cga awl 
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zeros(rawl,caw2); 

zeros(raw2,cag) zeros(raw2,cawl) aw2J; 

% 

%-B matrix : 

% 

[rbg,cbg] = size(bg}; 


[rbwl,cbwl] = size(bwl); 

[rbw2,cbw2] = size(bw2); 

bl = [zeros(rbg,cbwl);bwl;zeros(rbw2,cbwl)]; 
b2 = [bg;-bwl*dga;bw2]; 

% 

%-C matrix : 

% 

[rcgb^ccgb] = size(cgb); 

[rcwl,ccwl] = size(cwl); 

[rcw2,ccw2] = size(cw2); 

cl = [-dwl*cga cwl zeros(rcwl,ccw2);zeros(rcw2,ccgb) 
zeros(rcw2,ccwl) cw2; 

cgb zeros(rcgb,ccwl) zeros(rcgb , ccw2)]; 

[rcga,ccga] = size(cga); 

c2 = [-cga zeros(rcga,ccwl) zeros(rcga,ccw2)]; 

% 

%-D matrix : 


% 

[rdgb,cdgb] = size(dgb); 
[rdwl,cdwl] = size(dwl); 
[rdw2,cdw2] = size(dw2); 


dll = [dwl;zeros(rdw2/cdwl);zeros(rdgb,cdwl)]; 
dl2 = [-dwl*dga;dw2;dgb]; 

[rdga,cdga] = size(dga); 
d21 = eye(rdga); 
d22 = -dga; 

% 

%-End of AUGPL.M- % 
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%£18analy.in 
diary fl8analy.dat 
format short e 
disp(' *) 

dispC •) 

disp(* This script file analyzes the results of Hinf synthesis 
of the') 

dispC F-18 MIMO model by plotting resultant system return 

difference*) 

dispC singular values, calculating the poles/zeros,and making 
Bode and') 

disp(' Nyquist plots.') 
disp(* ') 

disp(* *) 

dispC The unaugmented state space plant representation*) 

dispC *) 

disp(' ') 

ag,bg,cg,dg 

disp(' ') 

dispC ') 

%pause 

dispC State space representation of the controller') 

disp(' ') 

dispC *) 

acp,bcp,ccp,dcp 

disp(' *) 

disp(* ') 

%pause 

dispC Open loop state space representation of controller') 

disp(' plant series') 

dispC *) 

dispC ') 

algf,blgf,clgf,dlgf 

[alfg,blfg,clfg,dlfg]=series(ag,bg,eg,dg,acp,bcp,ccp,dcp) 
disp(' ') 

disp(' *) 

dispC Compute singular value plot of return difference 

matrices') 
dispC ') 

disp(* ') 

w=logspace(-3,3,100); 

svg=sigma(ag,bg,cg,dg,3,w); svg=20*logl0(svg); 
semilogx(w,svg) 

title('F-18 Limited Performance SV PLOT (I + G)') 
xlabelCFREQUENCY - rad/sec') 
ylabel(*SV - dB') 
grid 

meta fl8svlL 
pause 

svgf=sigma(algf,blgf,clgf,dlgf,3,w); svgf=20*logl0(svgf); 
semilogx(w,svgf) 

title('F-18 Limited Performance SV PLOT (I + GF)') 







xlabel (•FREQUENCY - rad/sec') 

ylabel C SV - dB') 

grid 

meta fl8sv2L 
pause 

svfg=sigina(alfg,blfg,clfg,dlfg,3,w); svfg=20*logl0(svfg); 
semllogx(w,svfg) 

titleCF-18 Limited Performance SV PLOT (I + FG)') 
xlabel(•FREQUENCY - rad/sec') 
ylabelCSV - dB') 
grid 

meta fl8sv3L 
pause 

svgf=sigma(algf,blgf,clgf,dlgf,4,w); svgf=20*logl0(svgf); 
semilogx(w,svgf) 

title('F-18 Limited Performance SV PLOT (I + inv(GF))') 
xlabel('FREQUENCY - rad/sec') 
ylabelCSV - dB') 
grid 

meta fl8sv4L 
pause 

Syfg=sigma(alfg,blfg,clfg,dlfg,4,w); svfg=20*logl0{svfg); 
sem.ilogx (w, svfg) 

title CF-18 Limited Performance SV PLOT (I + inv(FG))') 
xlabel('FREQUENCY - rad/sec') 
ylabelCSV - dB') 
grid 

meta fl8sv5L 
pause 

svloop=sigma(algf,blgf,clgf,dlgf,l,w); svloop=20*logl0(svloop); 
semilogx(w,svloop) 

title('F-i8 Limited Performance SV PLOT (GF)') 
xlabel('FREQUENCY - rad/sec') 
ylabelCSV - dB') 
grid 

meta fl8sv6L 
pause 

dispC ') 

dispC ') 

disp(' Closed loop state space representation of 

controller/plant series *) 
dispC ') 

disp(' ') 

acgf , bcgf,ccgf,dcgf 
tzerocgf=t 2 ero(acgf , bcgf , ccgf,dcgf) 
dispC ') 

disp C ') 

%pause 

dispC Poles and zeros of the closed loop controller/plant 

system ') 

dispC (acgf, bcgf, ccgf, dgcf) ') 

dispC ') 
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dispC ') 

disp(* output alpha/g vs input 1 ') 
disp(* ') 

[z,p,k]=ss2 zp(acgf,bcgf , ccgf,dcgf ,1) 

%pause 
dispC •) 

disp(' output alpha/q vs input 2 ') 
dispC ') 

[z,p,k]=ss2 zp(acgf,bcgf , ccgf , dcgf ,2 ) 
pause 

dispC ') 

dispC ') 

disp(' Open loop Bode plots of outputs vs inputs') 

dispC •) 

disp C ') 

[inagl,phasel]=bode(algf ,blgf ,clgf ,dlgf,l,w); inagl=20*logl0(inagl}; 
[mag2,phase2]=bode(algf ,blgf ,clgf ,dlgf ,2,w) ; inag2=20*logl0(iaag2) ; 
seinilogx(w,inagl (:, 1)) 

title (’F-18 LIMITED PERFORMNCE OPEN LOOP BODE PLOT INPUT 1 / alpha 
(GF)') 

xlabelCFREQUENCY - rad/sec') 

ylabelCGAIN - dB') 

grid 

neta flSbodelL 
pause 

seinilogx(w,inagl(: ,2)) 

title C F-18 LIMITED PERFORMANCE OPEN LOOP BODE PLOT INPUT 1 / q 
(GF)*) 

xlabelCFREQUENCY - rad/sec') 

ylabelCGAIN - dB') 

grid 

neta fl8bode2L 
pause 

senil6gx(w,nag2(:,1)) 

title C F-18 LIMITED PERFORMANCE OPEN LOOP BODE PLOT INPUT 2 / alpha 
(GF)') 

xlabel('FREQUENCY - rad/sec') 

ylabelCGAIN - dB') 

grid 

neta fl8bode3L 
pause 

senilogx(w,nag2(: ,2)) 

titleCF-18 LIMITED PERFORMANCE OPEN LOOP BODE PLOT INPUT 2 / q 
(GF)') 

xlabel('FREQUENCY - rad/sec') 

ylabelCGAIN - dB') 

grid 

neta f jL8bode4L 
dispC ') 

dispC ') 

dispC Closed loop Bode plots of outputs vs inputs') 

dispC ') 
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disp(* ') 

[inagl,phasel]=bode(acgf ,bcgf ,ccgf ,dcgf,l,w); inagl=20*logl0(magi); 
[mag2,phase2]=bode(acgf ,bcgf ,ccgf ,dcgf ,2,w); mag2=20*]Logl0(mag2); 
semilogx(w,magl(:,!)) 

title('F-IS LIMITED PERFORMANCE CLOSED LOOP BODE PLOT INPUT 1 / 
alpha ') 

xlabel ('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

grid 

meta flSbodeSL 
pause 

semilogx(w,magl(:,2)) 

title('F-18 LIMITED PERFORMANCE CLOSED LOOP BODE PLOT INPUT 1 / q 

') 

Xlabel ('FREQUENCY - rad/sec') 

ylabelC'GAIN - dB') 

grid 

meta fl8bode6L 
pause 

semilogx(w,mag2(:,1)) 

title('F-18 LIMITED PERFORMANCE CLOSED LOOP BODE PLOT INPUT 2 / 
alpha ') 

xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

grid 

meta fl8bode7L 
pause 

semilogx(w,mag2(:,2)) 

title('F-18 LIMITED PERFORMANCE CLOSED LOOP BODE PLOT INPUT 2 / q 

') 

Xlabel ('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

grid 

meta fl8bode8L 
%pause 
disp(' ') 

%disp(' ') 

%disp(' Open loop nyquist plots of outputs vs inputs') 

%disp(' ') 

%disp(' ') 

%w=logspace(-4,5/200); 

%[rel,iml]=nyquist(algf,blgf,cigf,dlgf,l,w); 

%[re2 , im2]=nyquist(algf , blgf , clgf , dlgf,2,w); 

%rel=[rel;rel]; iml=[iml;-iml]; 

%re 2 =[re 2 ;re 23 ; im2=[im2;-im2]; 

%plot(rel(;,1),iml(:,1)) 

%title('NYQUIST PLOT INPUT 1 / alpha (GF)') 

%xlabel('real (GF)') 

%ylabel (' imag (GF)') 

%grid 

%meta x29nyqu 
%pause 
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%plot(rel(: ,2) ,2)) 

%title('NYQUIST PLOT INPUT 1 / q (GF)') 
%xlabel('real <GF)') 

%ylabel('iinag (GF)') 

%grid 

%meta 

%pause 

%plot(re2(; ,,1) ,iin2 ( :, 1)) 

%title('NYQUIST PLOT INPUT 2 / alpha (GF)•) 
%xlabel('real (GF)‘) 

%ylabel('imag (GF)') 

%grid 

%ineta 

%pause 

%plot (re2 {i ,2) I iin2 (:, 2)) 

%title('NYQUIST PLOT INPUT 2 / q (GF)') 
%xlabel('real (GF)') 

%ylabel('imag (GF)') 

%grid 

%meta 

%pause 

end 
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%fl8resp.in 
diary flSresp.dat 
format short e 
disp(' ') 

disp(' ') 

disp(' This script file is designed to calculate and plot the 

time•) 

dispC response of the augmented F-18 (controller/plant closed 
loop series)*) 

disp(* to 1 degree / 1 second step input from each of *) 
disp(' the two inputs separately.') 
disp(' ') 

disp(* ') 

disp(' Setting up the time vectors') 
disp(' ') 

time=[0:0.01:4]; 

%stas=[l 000 0;0 100 0;0 010 0;0 001 0;0 0001] 
%ccgfwc=[zeros(5,18) stas zeros(5,6)] 

%dcgfwc=zeros(5,2) 

%stas=[l 000000 0;0 100000 0;0 010000 0;0 00100 

0 0;... 

%0 0 0 0 1 0 0 0;0 000010 0;0 000001 0;0 0000001] 
%ccgfwc=- [ zeros (8,15) stas zeros (8,6)] 

%dcgfwc=zeros(8,2) 
disp(* ') 

disp(' Plotting alpha, q, ds, dlf, and dtf responses to a 

0.01745 rad / 1 ') 

disp(' second Step from input one') 

disp(' *) 

u=[0.01745*ones(l,101) zeros(l,300);zeros(l,401)]*; 

%[y]=lsim(acgf,bcgf,ccgfwc,dcgfwc,u,time) 

[y/X]=lsim(acgf,bcgf,ccgf,dcgf,u,time); 

%plot(time,y(:,1)) 

%title('F-18 alpha RESPONSE TO 6.01745 rad / 1 sec STEP (INPUT 1) 
W2=-.001*) 

%xlabel('TIME - sec') 

%ylabel('radians') 

%grid 

%meta fl8rspllU 
%pause 

%plot(time,y(:,2)) 

%title('F-18 q RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 1) 
W2=-.001') 

%xlabel('TIME - sec') 

%ylabel('radians/second') 

%grid 

%meta fl8rspl2U 
%pause 

plot(time,x(:,17)) 

title('F-18 DS FOR 0.01745 rad / 1 sec STEP (INPUT 1) W2=-.001') 
xlabel('TIME sec') 
ylabel('radians') 




grid 

%ineta fl8rspl3L 
pause 

plot(time,y(:,1),timely(:,2)) 

title('F18 OPT. RESPONSE TO 1 DEG / ONE SECOND IMPULSE (INPUT 1) ') 

xlabelCTIME -SEC) 

ylabel('RADIANS') 

text(1.2,.0102,'ALPHA RESPONSE') 

text(0.3,.0012,'q RESPONSE') 

grid 

meta alphaqlT2U 
pause 

plot(tiTne,x(:,21)) 

title('F-18 DLF FOR 0.01745 rad / 1 sec STEP (INPUT 1) W2=-.001') 
XlabelCTIME - sec') 
ylabel('radians') 
grid 

%ineta fl8rspl7L 
pause 

plot(time,x(:,23)) 

title('F-18 DTF FOR 0.01745 rad / 1 sec STEP (INPUT 1) W2=-.001') 
xlabel('TIME - sec') 
ylabel('radians') 
grid 

%meta fl8rspl9L 
pause 

plot(time,X(:,18)) 

title('F-18 DSDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1) W2=-.001') 
XlabelCTIME - sec') 
ylabel('radians/second') 
grid 

%meta fl8rspl4L 
pause 

plot(time,x(:,22)) 

title('F-18 Lim. DLFDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1) 
W2=-.001') 

xlabel('TIME - sec') 
ylabel('radians/second') 
grid 

%meta fl8rspl8L 
pause 

plot(time,x(:,24)) 

title('F-18 Lim. DTFDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1) 
W2=-.OOr') 

Xlabel('TIME - sec') 
ylabel('radians/second') 
grid 

%meta fl8rspllOL 
pause 

disp(' ') 

dispC Plotting alpha, q, dc, df, and ds responses to a 0.01745 
rad / 1') 
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disp(• second step input from input two') 
disp(' ') 

u=[zeros(l,401);0.01745*ones(1,101) zeros(1,300)]*; 

% ty]=lsim(acgf,bcgf,ccgfwc,dcgfwc , u,time) 

[y£» X] =lsim (acgf, bcgf, ccgf, dcgf, u, time) 7 
plot(time,y(: , 1 ) ,time,y(: , 2 )) 

title('F-IS OPT. RESPONSE TO 0.01745 rad/ 1 sec STEP(INPUT 2)') 

xlabel('TIME - sec') 

ylabel('radians') 

text(1.2,.011,'q RESPONSE') 

text(1.2,-.001,'ALPHA RESPONSE') 

grid 

meta alphaq212U 
pause 

%plot(time,y(:, 2 )) 

%title('F-18 q RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 2) 
W2=-.001') 

%xlabel('TIME - sec') 

%ylabel('radians/second') 

%grid 

%meta fl8rsp22U 
%pause 

plot(time,x(:,17)) 

title('F-18 DS FOR 0.01745 rad / 1 sec STEP (INPUT 2 ) W2=-.001') 
xlabelCTIME - sec') 
ylabel('radians') 
grid 

%meta fl8rsp23L 
pause 

plot(time,x(:, 21 )) 

title('F-18 DLF FOR 0.01745 rad / 1 sec STEP (INPUT 2) W2=-.00l') 
xlabel('TIME - sec') 
ylabel('radians') 
grid 

%meta fl8rsp27L 
pause 

plot(time,x(;,23)) 

title('F-18 DTF FOR 0.01745 rad / 1 sec STEP (INPUT 2) W2=-.001') 
XlabelCTIME - sec') 
ylabel('radians') 
grid 

%meta fl8rsp29L 
pause 

pl6t(time,x(:,18)) 

title('F-18 DSDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2) W2=-.001') 
'•label (' TIME - sec') 

7 label('radians/second') 
c rid 

%meta fl8rsp26L 
pause 

plot(time,x(;, 22 )) 

title ('F-18 DLFDOT FOR 0.01745 rad / 1 seC STEP (INPUT 2) 
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W2=-.001') 

xlabelCTIME - sec') 
ylabeT('radians/second') 
grid 

%meta fl8rsp27L 
pause 

plot (tiine,x(: ,24)) 

title ('F-18 DTFDOT FOR 0.01745 rad / 1 sec STEP 

W2=-.001') 

xlabel('TIME - sec') 
ylabel('radians/second') 
grid 

%iaeta fl8rsp210L 
end 


(INPUT 2) 


94 



APPENDIX B 


F-18 STATE SPACE MODELS 

1) Open Loop F-18 Plant 

2) Optimum Hoo Closed Loop Controller 

3) Limited Performance H® closed Loop Controller 
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% Open Loop Plant of the F-18 
>> ag 


ag = 

Columns 1 through 5 
-1.1400e+00 9.7938e-01 

-7.0738e+00 -8.0902e-01 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


-9.9461e-02 

-7.9063e+00 

0 

0 

0 

-1.4691e+03 

0 

0 

0 

0 


0 

0 

l.OOOOe+00 

0 

0 

-4.9559e+01 

0 

0 

0 

0 


0 

0 

0 

l.OOOOe+00 

0 

-1.6122e+00 

0 

0 

0 

0 


Columns 6 through 10 


0 

0 

0 

0 

l.OOOOe+04 

-1.5410e+02 

0 

0 

0 

0 


5.6206e-04 

-1.7235e+00 

0 

0 

0 

0 

0 

-2.2301e+03 

0 

0 


0 

0 

0 

0 

0 

0 

l.OOOOe+00 

-1.0980e+02 

0 


3.9669e-03 
1.8786e+00 
0 
0 
0 
0 
0 
0 
0 

0 -1.2250e+03 


0 

0 

0 

0 

0 

0 

0 

0 

l.OOOOe+00 

-4.9700e+01 


» bg 
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bg = 


0 
0 
0 

2.1377e+03 
-3.0532e+05 
2.7277e+03 
0 
0 
0 
0 

» eg 
eg = 

1 0 0 0 0 0 0 

0 1 0 0 0 0 0 

» dg 
dg = 

0 0 

0 0 

» % Poles of the F-18 Open Loop Plant 
» eig(ag) 


0 

0 

0 

0 

0 

0 

0 

2.2301e+03 

0 

1.2250e+03 


ans = 


^9.7451e-01 

-9.7451e-01 

-6.2126e+01 

-6.2126e+01 

-1.4924e+01 

-1.4924e+01 

-2.6902e+01 

-8.2898e+01 

-2.4850e+01 

-2.4850e+01 


+ 2.6269e+00i 

- 2.6269e+00i 
+ 8.5022e+0li 

- 8.5022e+01i 
+ 3.3i99e+0li 

- 3.3i99e+0li 


+ 2.4647e+0li 
- 2.4647e+01i 


» diary off 


0 

0 


0 

0 
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» % Optimal Closed Loop State Space F-18 H-Infinity Controller 
» acgf 

acgf = 


Columns 1 through 6 


-8.5374e+04 

-1.0132e+04 

-3.1614e+03 

1.4349e+04 

-6.2351e+00 

3.6195e+03 

-1.0476e+05 

-1.2541e+04 

-3.8866e+03 

1.7514e+04 

-7.5004e+00 

4.2957e+03 

2.2082e+06 

2.6104e+05 

8.1683e+04 

-3.7082e+05 

1.6204e+02 

-9.5204e+04 

2.0282e+02 

-8.9659e+02 

-3.0481e+01 

-2.0145e+02 

2.4393e-01 

3.6679e+02 

-1.1859e-01 

-6.6984e-01 

-4.0737e-02 

-1.9136e-01 

-6.0713e-04 

-2.6865e-01 

-5.9038e+04 

-7.7319e+03 

-2.2235e+03 

9.7030e+03 

-4.0591e+0p 

2.3185e+63 

2.2799e+02 

1.3562e+03 

8.0283e+01 

3.7452e+02 

-7.9406e-01 

4.2695e+02 

-1.1229e+00 

-6.9592e+00 

•^4.0987e-01 

-1.9229e+00 

4.0321e-03 

-2.1645e+00 

-1.6276e+04 

-2.5326e+03 

-6.3355e+02 

2.5765e+03 

-9.5432e-01 

5.6803e+02 

4.5355e-02 

3.8317e-03 

1.5852e-03 

-8.4079e-03 

2.2515e-05 

-3.2101e-03 

-2.6327e-03 

2.0579e-03 

4.8972e-05 

1.5237e-03 

-8.1576e-06 

2.1350e-03 

5.1076e-03 

4,2810e-03 

4.1494e—04 

7.3066e--04 

8i0206e-06 

2.8272e-03 

1.0102e-04 

2.2190e-04 

1.7064e-05 

9.3235e-05 

-1.7261e-06 

1.9550e-04 

-5.3128e-03 

-1.3122e-02 

-9.8395e-04 

-5.1009e-03 

-3.3955e-06 

-1.1396e-02 

0 

n 

0 

0 

0 

0 

0 

n 

0 

0 

0 

0 

u 

0 

0 

0 

0 

0 

u 

-1.5489e+04 

-1.8299e+03 

^5.7521e+02 

2.6270e+03 

-1.1398e+00 

6.1035e+02 

2.2122e+06 

2.6135e+05 

8.2156e+04 

-3.7521e+05 

1.6280e+02 

-8.7174e+04 

-1.9764e+04 

-2.3349e+03 

-7.3397e+02 

3.3521e+03 

-1.4544e+00 

7.7880e+62 

0 

A 

0 

0 

0 

0 

U 

-5.9035e+03 

-5.8679e+02 

-2.1394e+02 

1.0604e+03 

4.1411e-02 

2.6132e+02 
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0 

0 

0 

0 

0 


~3.2428e+03 -3.2233e+02 

l._4355e+02 

Columns 7 through 12 

-1.17526+02 

5.82486+02 

2.27476-02 

* 

1.6525e+04 

-7.4498e+05 

-8.12136+01 

-2.32316+04 

-8.53926+05 

6.92076+06 


2.0288e+04 

-1.4361e+06 

-9.97006+01 

-2.80966+04 

-1.03506+06 

8.36156+06 


“4.2803e+05 

1.7108e+07 

2.10366+03 

6.00496+05 

2.21326+07 

-1.79506+08 


-3.1957e+02 

2.3416e+05 

1.73426+00 

-7.65476+02 

-5.87526+03 

5.63016+04 


-1.2503e-01 

-9.4379e+02 

6.82866-04 

3.75896-01 

2.09666+01 

-2.18366+02 


1.1240e+04 

-8.0642e+05 

-5.51796+01 

-1.58216+04 

-5.82706+05 

4.71626+06 


2.59156+02 

1.85206+06 

-1.43096+00 

-5.71136+02 

-4.04386+04 

4.26106+05 


-1.34906+00 

-9.31526+03 

6.43426-03 

2.87406+00 

2.03746+02 

-2.14446+03 


2.98826+03 

-4.57096+05 

-1.46256+01 

-4.34136+03 

-1.55286+05 

1.24596+06 

1 # 

-9.50546-03 

-1.93586+01 

3.96276-05 

1.47796-02 

-1.08946+01 

6.41796+01 


1.16396-03 

-2.66266+00 

-3,28696-06 

-4.43386-03 

-3.71776+00 

-2-58026+01 

- 

-7.73866-06 

-2.19526+01 

8.63956-07 

-4.06366-03 

1.18246+00 

2.43806+00 


4.35516-05 

8.44196-02 

-4.54446-08 

-3.55376-04 

-1.78686+00 

-3.34366+00 


-2.36566-03 

1.80886+01 

9.69506-06 

2-00696-02 

-7.27026-01 

2.81986+00 


0 

n 

0 

0 

0 

0 


0 

n 

0 

0 

0 

0 


0 

n 

0 

0 

0 

0 


w 

2.98586+03 

-1.20316+05 

-1.46836+01 

-4.23826+03 

-1.55296+05 

1.25956+06 


-4.26456+05 

1.71846+07 

2.09716+03 

6.05336+05 

2.21806+07 

-1.79896+08 

> 

3.80996+03 

-1.53526+05 

-1-87356+01 

-5.40796+03 

-1.98156+05 

1.60716+06 


0 

n 

0 

0 

0 

0 

• 

1.18736+03 

-1.82306+06 

-5.75976+00 

-1.70016+03 

-1.93066+04 

6.81546+04 


0 

0 

0 

0 

0 
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0 

6.5219e+02 

-3.1638e+00 

-9.33846+02 

-1.06056+04 

3.74376+04 

1.0014e+06 

Columns 13 through 18 

7.0693e+07 -8.8082e+07 

r\ 

6.38516-06 

5.61216-07 

0 

8.6707e+07 

-9.2874e+07 

1.74896-05 

1.22536-06 

0 

u 

-1.8283e+09 

2.J406e+09 

1.27406-06 

9.42376-08 

0 

u 

-7.1948e+03 

-6.7861e+06 

6.49876-06 

5.09626-07 

0 

u 

1.0692e+02 

2.7226e+04 

2.14916-08 

-2.87176-09 

0 

u 

4.8902e+07 

n 

-5.2479e+07 

1.42236-05 

1.07176-06 

0 

u 

-2.3186e+05 

A 

-5.3390e+07 

3.28736-06 

2.40956-07 

0 

u 

1.1552e+03 

n 

2.68566+05 

1.77946-08 

-7.84336-09 

0 

U 

1.3509e+07 

n 

-7.6963e+06 

-2.46796-05 

-1.93796-06 

0 

u 

6.7517e+02 

n 

-5.6786e+02 

1.53316-03 

-1.07276-02 

0 

U 

-2.7816e+02 

n 

3-01496+02 

-1.46566-02 

-4.87336-02 

0 

u 

7.1802e+01 

n 

5.24046+02 

-6.80266-02 

1.80076-02 

0 

u 

-4.1241e+01 

2.80426-02 

1.38206-01 

5.55726+00 

0 

5.9266e-02 

-4.66536+02 

-2.12996+01 

5.19506-01 

0 

u 

0 

0 

-1.14006+00 

9.79386-01 

-9.94616-02 

U 

0 

n 

0 

-7.07386+00 

-8.09026-01 

-7.90636+00 

u 

0 

0 

0 

0 

0 

l.OOOOe+00 

1.2829e+07 

n 

-1.64166+07 

0 

0 

0 

u 

-1.8323e+09 

A 

2-34466+09 

0 

0 

0 

W 

1.6369e+07 

-2.09466+07 

0 

0 

-1.46916+03 

4.9559e+01 

0 

A 

0 

0 

0 

0 

U 

5.0852e+06 

A 

4.50206+07 

0 

0 

0 

U 

0 

0 

0 

0 

0 

0 
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2.7933e+06 2.4730e+07 

0 

Columns 19 through 24 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5.6206e-04 

0 

-3.9669e-03 

0 

0 

0 

-1.7235e+00 

0 

-1.8786e+00 

0 

0 

0 

0 

0 

0 

l.OOOOe+00 

0 

0 

0 

0 

0 

0 

0 

l.OOOOe+04 

0 

0 

0 

-1.6122e+00 

0 

-i.5410e+02 

0 

0 

0 

0 

0 

0 

0 

l.OOOOe+00 

0 

0 

0 

0 

-2.2301e+03 

-1.0980e+02 

0 

0 

l.OOOOe+00 

0 

0 

0 

0 

0 

0 

0 

0 

-1.2250e+03 



-4.9700e+01 


» bcgf 


bcgf = 


-6.3851e-06 

-5.61216-07 

-1.74896-05 

-1.22536-06 

-1.27406-06 

-9.42376-08 

-6.49876-06 

-5.09626-07 

-2.14916-08 

2.87176-09 

-1.42236-05 

-1.07176-06 

-3.28736-06 

-2.40956-07 

-1.77946-08 

7.84336-09 

2.46796-05 

1.93796-06 

-1.53316-03 

1.07276-02 

1.46566-02 

4.87336-02 

6.80266-02 

-1.80076-02 

-1.38206-01 

-5.55726+00 

2.12996+01 

-5.19506-01 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

» ccgf 


ccgf = 


Columns 1 through 12 

0 0 

0 0 

0 


0 0 

0 0 

0 


Columns 13 

through 24 

0 0 

1 0 

0 



0 0 0 
0 0 0 

0 0 0 


0 0 0 
0 0 0 

0 0 0 


c 


t 


102 





» % Limited Performance F-18 H-Infinity Closed Loop Controller 
» acgf 

acgf = 

Columns 1 through 6 


5.4853e+01 

2.7764e+00 

6.8308e+01 

-2.4157e~03 

—4.5556e+01 

-1.3317e-03 

-8.2591e+02 

-1.2620e+01 

-7.1443e+02 

-1.2218e-01 

4.5375e+02 

3.9671e-02 

5.1971e-01 

3.2926e-03 

3.5172e-01 

-1.7881e-02 

-2.4551e-01 

-3.0125e-05 

-9.5636e+02 

-5.8066e+00 

-6.4519e+02 

-2.2369e-01 

4.5000e+02 

5.6299e-02 

-2.3815e-01 

-1.4021e-03 

-1.5991e-01 

-5.5533e-05 

1.1168e-01 

-1.7986e-02 

9.6310e+01 

3.1452e+02 

8.8451e+03 

-3.9933e+00 

-2.7804e+03 

5.3889e-01 

9.9182e+01 

2.3513e+00 

1.0866e+02 

-1.1930e-02 

-6.1211e+01 

5.3155e-04 

9.5922e+01 

-2.9221e+00 

-2.0037e+01 

7.2362e-02 

-1.7863e+01 

-1.4373e-02 

-1.5046e+00 

-1.6638e-01 

-3.3465e+00 

1.3703e-03 

2.3553e+00 

-3.2388e-04 

2.6692e-02 

-9.7550e-03 

-1.7103e-01 

1.9213e-04 

1.1688e-01 

-4.0487e-05 

3.1661e-03 

-1.3737e-04 

-1.2917e-03 

-1.4128e-06 

-1.1759e-03 

2.4284e-07 

-5.01380-04 

1.2726e-05 

-2.2227e-04 

-1.8603e-06 

6.3267e-04 

6.1350e-07 

-1.6750e-03 

9.4405e-07 

-9.6082e-04 

-1.4759e-06 

1.4166e-03 

8.2410e-07 

-9.0621e-03 

3.7032e-05 

-4.9989e-03 

-5.5652e-06 

7.9886e-03 

3.9378e-06 

0 

n 

0 

0 

0 

0 

u 

0 

n 

0 

0 

0 

0 

u 

0 

n 

0 

0 

0 

0 

u 

1.4537e+00 

2.9299e+O0 

5.2614e+01 

-3i4809e-02 

-1.4738e+01 

6.7451e-03 

-2.0763e+02 

-4.1847e+02 

-7.5146e+03 

4.9716e+00 

2.1050e+03 

-9.6337e-01 

1.8550e+00 

3.7386e-i00 

6.7135e+01 

-4.4416e-02 

-1.8805e+01 

8.6067e-03 

0 

n 

0 

0 

0 

0 

u 

8.2662e+00 

1.3263e+00 

2.5557e+01 

-5.3459e-03 

-1.1852e+01 

1.3675e-03 

0 

0 

0 

0 

0 
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0 

4.5406e+00 

1.40396+01 

-2.93656-03 

-6.51026+00 

7.51166-04 

> 

7.2854e-01 

Columns 7 through 12 

2.42786+01 2.19706+02 

-7.44806+02 

3.45316+02 

-1.28836+03 

A 

-3.37646+03 

1.14756+02 

-1.52756+03 

-7.44946+02 

-1.20816+03 

4.42616+04 


—9.22656+04 
-9.51756-02 

-1.79716-01 

5.67856-01 

-2.90316-01 

1.10636+00 


2.09116+00 

1.69486+02 

3.55406+02 

-1.21436+03 

4.96146+02 

-1.68246+03 


-5.27226+03 

4.42746-02 

9.27006-02 

-2.99616-01 

1.47586-01 

-5.43026-01 


-1.18086+00 

-1.91256+03 

1.15616+05 

6.22906+04 

7.51736+04 

-3.36426+06 


7.19216+06 

-5.84296+01 

5.43036+02 

-2.72876+02 

-3.38316+02 

-1.24676+04 


2.78436+04 

5.84736+01 

-1.22806+03 

-2.70126+02 

-9.37986+02 

3.35466+04 


-6.87526+04 

1.18176+01 

-3.56956+01 

-1.14116+01 

-9.24596+01 

1.03206+03 


-1.72966+03 

6.99026-01 

-2.72566+00 

2.75616+00 

-2.50036+00 

7.51736+01 


-1.29756+02 

-2.10936-03 

-5.02106-02 

-5.96446-02 

-5.48256-02 

-1.89136+01 


4.24886+01 

1.75456-03 

1.97976-03 

-2.08386-02 

1.54716-01 

-3.45656+00 


-8.83316+00 

3.35866-03 

2.93266-03 

-2.95696-02 

4.03746-01 

-1.17786+00 


-4.02336+00 

1.88706-02 

1.90256-02 

-1.71116-01 

2.11686+00 

1.16946+00 


-5.78766-02 

0 

0 

0 

n 

0 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


\J 

-8.22036+01 

8.04396+02 

4.30556+02 

5.26496+02 

-2.35566+04 


5.03596+04 

1.17416+04 

-1.14896+05 

-6.14946+04 

-7.51976+04 

3.36456+06 


-7.19266+06 

-1.04896+02 

1.02646+03 

5.49386+02 

6.71806+02 

-3.00586+04 

i 

6.42586+04 

0 

n 

0 

0 

0 

0 


-5.89086+01 

3.11326+02 

4.22186+02 

6.31366+02 

-1.06516+04 

>1 

1.37206+04 

0 

0 

0 

0 

0 

0 
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-3.2358e+01 

1.71016+02 

2.31916+02 

3.46816+02 

-5.85096+03 

7.5365e+03 

Columns 13 through 18 

1.5165e+04 3.6975e+05 

r\ 

9.17836-04 

-1.30666-03 

0 

U 

-3.9813e+05 

3.29316+05 

5.16166-04 

-7.33206-04 

0 

U 

-1.1531e+01 

r\ 

-2.55446+02 

4.90556-07 

-6.51766-07 

0 

U 

1.9814e+04 

5.31196+05 

-7.90866-04 

1.14346-03 

0 

u 

5.8190e+00 

1.35816+02 

-3.69446-07 

5.44846-07 

0 

u 

3.0288e+07 

n 

-3.47456+07 

-2.48936-06 

3.75706-06 

0 

U 

1.23530+05 

r\ 

-6.74876+04 

-1.83756-03 

2.72346-03 

0 

U 

-3.02930+05 

1.70516+05 

-1.76206-03 

2.56936-03 

0 

U 

-8.60930+03 

-8.27876+03 

1.63186-02 

-2.40676-02 

0 

u 

-6.65220+02 

A 

-6.40546+02 

-2.08636-01 

3.11586-01 

0 

u 

1.89300+02 

n 

-2.50116+01 

-1.53736-01 

2.41146-01 

0 

U 

-2.66890+01 

0 

-1.37930+01 

1.42046+01 

-4.27176-01 

2.57816-02 

0 

9.22626-01 

3.75636-01 

8.97726-01 

0 

U 

-5.9566e-01 

0 

0 

-2.14686+01 

-3.21206+00 

1.09466+00 

0 

0 

-1.14006+00 

9.79386-01 

-9.94616-02 

0 

0 

0 

-7.07386+00 

-8.09026-01 

-7.90636+00 

0 

0 

0 

0 

0 

0 

l.OOOOe+00 

2.1208e+05 

n 

-2.43286+05 

0 

0 

0 

U 

-3.0290e+07 

3.47476+07 

0 

0 

0 

U 

2.70616+05 

-3.10426+05 

0 

0 

-1.46916+03 

-4.95596+01 

0 

r\ 

0 

0 

0 

0 

U 

1.05526+05 

5.04926+05 

0 

0 

0 

u 

0 

0 

0 

0 

0 

U 

5.79606+04 

2.77356+05 

0 

0 

0 
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Columns 19 through 24 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5.6206e-04 

0 

-3.9669e-03 

0 

0 

0 

-i.7235e+00 

0 

-1.8786e+00 

0 

0 

0 

0 

0 

0 

l.OOOOe+00 

0 

0 

0 

0 

0 

0 

0 

l.OOOOe+04 

0 

0 

0 

-1.6122e+00 

0 

-1.5410e+02 

0 

0 

0 

0 

0 

0 

0 

l.OOOOe+00 

0 

0 

0 

0 

-2.2301e+03 

-1.0980e+02 

C 

0 

l.OOOOe+00 

0 

0 

0 

0 

0 

-4.9700e+01 

0 

0 

0 

-1.2250e+03 





>> bcgf 
bcgf = 


-9.17836-04 

1.30666-03 

-5.16166-04 

7.33206-04 

-4.90556-07 

6.51766-07 

7.90866-04 

-1.14346-03 

3.69446-07 

-5.44846-07 

2.48936-06 

-3.75706-03 

1.83756-03 

-2.72346-03 

1.76206-03 

-2.56936-03 

-1.63186-02 

2.40676-02 

2.08636-01 

-3.11586-01 

1.53736-01 

-2.41146-01 

4.27176-01 

-2.57816-02 

-3.75636-01 

-8.97726-01 

3.21206+00 

-1.09466+00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


» ccgf 
ccgf = 

Columns l through 12 

00000000000 

0 

00000000000 

0 


Columns 13 through 24 
0 0 10 

0 

0 0 0 1 

0 



\ 






108 




» dcgf 

dcgf = 

0 0 

0 0 

» % Poles of the F-18 Limited Performance Closed Loop Controller 
» eig(acgf) 

ans = 


-1.0017e+03 

-6.2126e+01 + 3,5022e-f0li 
-6.2126e+01 - 8.5022e+0li 
~1.7711e+01 + 8.1027e+01i 
-1.7711e+01 - 8.1027e+0li 
-8.2898e-f01 
-7.6597e+01 

-1.4924e+01 + 3.3199e+0li 
-1.4924e+01 - 3.3199e+0li 
-2.4200e+01 + 2.5750e+0li 
-2.42C0e+01 - 2.5758e+0li 
-2.4850e+01 + 2.4647e+0li 
-2.4850e+01 - 2.4647e+Cli 
-2.0234e+01 + 2.0813e+0li 
-2.0234e+01 - 2.0813e+0li 
-2.6902e+01 

-2.5668e+01 + 7.8886e-0li 
-2.5668e+01 - 7.8886e-0li 
-9.7451e-01 + 2.6269e+00i 
-9.7451e-01 - 2.6269e+Q0i 
-2.5001e+00 
-1.9977e+00 
-1.8000e-02 
-1.8000e-02 

» diary off 
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